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ABSTRAGT 

We present ALMA cycle 0 observations of the molecular gas and dust in the IR-bright mid-stage 
merger VV 114 obtained at 160 - 800 pc resolution. The main aim of this study is to investigate 
the distribution and kinematics of the cold/warm gas and to quantify the spatial variation of the 
excitation conditions across the two merging disks. The data contain 10 molecular lines, including the 
first detection of extranuclear GH3OH emission in interacting galaxies, as well as continuum emission. 

We map the (3-2)/^^CO (1-0) and the (1-0)/^^CO (1-0) line ratio at 800 pc resolution 
(in the units of K km s“^), and find that these ratios vary from 0.2 - 0.8 and 5 - 50, respectively. 
Conversely, the 200 pc resolution HCN (4-3)/HCO+ (4-3) line ratio shows low values (< 0.5) at a 
filament across the disks except for the unresolved eastern nucleus which is three times higher (1.34 
± 0.09). We conclude from our observations and a radiative transfer analysis that the molecular gas 
in the VV 114 system consists of five components with different physical and chemical conditions; 
i.e., 1) dust-enshrouded nuclear starbursts and/or AGN, 2) wide-spread star forming dense gas, 3) 
merger-induced shocked gas, 4) quiescent tenuous gas arms without star formation, 5) H2 gas mass 
of (3.8 ± 0.7) X 10^ Mq (assuming a conversion factor of aco = 0.8 Mq (K km s“^ pc^)“^) at the 
tip of the southern tidal arm, as a potential site of tidal dwarf galaxy formation. 

Subject headings: galaxies: individual (VV 114, IC 1623, Arp 236) — galaxies: interactions — galaxies: 
starburst — galaxies: nuclei — ISM: molecules 


1. INTRODUCTION 

Galaxy interactions and mergers play important roles 
in triggering star formation and/or fueling the nuclear 
activity in the merging host galaxies (Hopkins et al. 
2006). Recent high resolution simulations of major merg¬ 
ers show that large scale tidal forces as well as small scale 
turbulence and stellar feedback can significantly influ¬ 
ence the distribution of gas, forming massive clumps of 
dense gas with Mhj = 10® - 10® Mq (e.g., Teyssier et 
al. 2010; Hopkins et al. 2013). These simulations also 
predict that the star formation not only increases as the 
galaxies first collide, but it also persists at a higher rate 
throughout the merger process, peaking at the final coa¬ 
lescence. 

(Ultra-)Luminous Infrared Galaxies (U/LIRGs; Soifer 
et al. 1987) at low redshifts are almost exclusively 
strongly interacting and merging systems (Kartaltepe et 
al. 2010), often found at the mid to final stages of the 
merger. The elevated level of infrared luminosity origi¬ 
nates from the reprocessed emission from the dust par¬ 
ticles surrounding the starburst or the Active Galactic 
Nuclei (AGNs), both of which are likely triggered by 
the tidal interaction. The highest gas surface densities 
(Eh 2 = 5.4 X 10^ - 1.4 X 10® Mq pc“^) and conse- 
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quently the highest star formation activities (Xsfr = ~ 
1000 Mq yr“^ kpc“^) are usually found near the com¬ 
pact nuclear region (e.g., Arp22C), NGC 6240; Downes 
& Solomon 1998; Engel et al. 2010; Wilson et al. 2014). 
Dense molecular gas {n ~ 10® - 10^ cm“®) in U/LIRGs 
directly shows nuclear gas distribution and kinematics 
(e.g., lono et al. 2004; Sakamoto et al. 2014). They are 
often surrounded by diffuse gas (n ^ 10^ - 10® cm“®) 
that may or may not be directly associated with star 
formation activities. 

It has been demonstrated that the HGN (4-3) and 
HCO’*' (4-3) emission lines, whose critical densities are 
8.5 X 10® and 1.8 x 10® cm“®, respectively, can be used 
as tracers of the dense gas (e.g., lono et al. 2013; Garcia- 
Burillo et al. 2014; Imanishi & Nakanishi 2014). On 
the other hand, GO (1-0) and ^®CO (1-0) line emis¬ 
sion, whose critical densities are 4.1 x 10^ and 1.5 x 10® 
cm“®, respectively, have been used extensively for trac¬ 
ing the global gas distribution and kinematics in merging 
U/LIRGs (e.g., Yun et al. 1994; lono et al. 2004; Ueda 
et al. 2014). In addition, the ratio of these lines (e.g. 
®^GO/®®CO and HGN/HGO+) have been used to inves¬ 
tigate the properties of the ISM (Gasoli et al. 1992; Aalto 
et al. 1997) or to search for buried AGNs (e.g., Imanishi 
et al. 2007; Imanishi & Nakanishi 2014). Limitations in 
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Fig. 1.— The HST/ACS image of VV 114 (Evans 2008). There 
is a dust lane from north to south in front of the eastern galaxy. 
The red crosses show the positions of the nuclei defined by the peak 
positions of the Ks-band observation (Tateuchi et al. 2012). The 
white ellipse shows a field of view of 3-point mosaic observation 
with band 7, while the white circle shows a field of view of 7-point 
mosaic observation with band 7 (see §2). 


sensitivity and angular resolution have been the major 
obstacles in understanding the detailed distribution and 
kinematics of both dense and diffuse gas, and investigat¬ 
ing the spatial variation of the line ratios and the physical 
condition of gas. 

In this paper, we present Atacama Large Millime¬ 
ter/submillimeter Array (ALMA) cycle 0 observations 
of the IR-bright merging galaxy VV 114. VV 114 is one 
of the best samples for studying the gas response during 
the critical stage when the two gas disks merge (lono et 
al. 2005; Wilson et al. 2008). The target molecular lines 
include (1-0), (1-0), (3-2), HCN (4- 

3) and HCO+ (4-3), and we also present the maps of 
CH 3 OH (2fc-lfc), CS (2-1), CN (I 1 / 2 -O 1 / 2 ), CN (I 3 / 2 - 
O 1 / 2 ), and CS (7-6) lines which were observed simulta¬ 
neously within the same band. The main aim of this 
study is to investigate the distribution and kinematics of 
the diffuse and dense molecular gas and to quantify the 
spatial variation of the excitation conditions across the 
two merging disks. 

VV 114 is a gas-rich (Mhj = 5.1 x 10^° Mq; Yun 
et al. 1994) nearby (D = 82 Mpc; I'.'O = 400 pc) inter¬ 
acting system (Figure 1) with high-infrared luminosity 
(Lir = 4.7 X 10^^ Lq; Armus et al. 2009). The pro¬ 
jected nuclear separation between the two optical galax¬ 
ies (VV 114E and VV 114W) is about 6 kpc. Frayer et 
al. (1999) found a large amount of dust (Mdust = 1-2 
X 10® Mq) distributed across the system with a dust 
temperature of 20 - 25 K. About half of the warmer 
dust traced in the mid-lR (MIR) is associated with the 
eastern galaxy, where both compact (nuclear region) and 
extended emission is found (Le Floc’h et al. 2002). Rich 
et al. ( 2011 ) found a bimodal distribution of velocity dis¬ 
persions of several atomic forbidden lines and emission 
line ratios indicative of composite activity explained by 
a combination of wide-spread shocks and star formation. 
The wide-spread star formation is also revealed by Paa 
observation using ANIR camera mounted on miniTAO 
(Tateuchi et al. 2012, see also Appendix A.l). lono et 


al. (2013) (hereafter paper 1) identified a highly obscured 
AGN and compact starburst clumps using sub-arcsecond 
resolution ALMA cycle 0 observations of HCN (4-3) and 
HCO'*' (4-3) emission. 

This paper is organized as follows. We describe our ob¬ 
servations and data reduction in §2, and results in §3. In 
§4 and §5, we provide molecular line ratios and physical 
parameters, such as the gas/dust mass, the gas tempera¬ 
ture, and the gas density. In § 6 , we present the properties 
of “dense” gas (§ 6 . 1 ), the comparison between molecular 
gas and star formation (§6.2), the discussions of the CO 
isotope enhancement (§6.3), the gas-to-dust mass ratio 
(§6.4), the fractional abundances of CS, CH 3 OH, and CN 
relative to H 2 (§6.5), and a potential tidal dwarf galaxy 
formations at the tip of the tidal arms of VV 114 (§ 6 . 6 ). 
We summarize and conclude this paper in §7. Through¬ 
out this paper, we adopt Hq = 73 km s“^ Mpc“^, Om = 
0.27, and Oa = 0.73. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations toward VV 114 were carried out as an 
ALMA cycle 0 program (ID = 2011.0.00467.S; PI = D. 
lono) using fourteen - twenty 12 m antennas. The band 3 
and band 7 receivers were tuned to the ^^CO (1-0), 
CO (1-0), i^CO (3-2), HCN (4-3), and HCO+ (4- 
3) line emissions in the upper side band (see Table 1). 
The ^^CO (1-0) data were obtained on November 6 , 2011 
and May 4, 2012 in the compact and extended configura¬ 
tions, respectively. The ^®CO (1-0) data were obtained 
on May 27 and July 2, 2012 in the compact configuration. 
The ^^CO (3-2) emission was observed on November 5, 
2011 in the compact configuration (7-point mosaic). The 
HCN (4-3) and HCO+ (4-3) data were obtained on July 
1, 2, and 3, 2012 in the extended configuration (3-point 
mosaic), simultaneously. Each spectral window had a 
bandwidth of 1.875 GHz with 3840 channels, and two 
spectral windows were set to each sideband to achieve 
a total frequency coverage of ^ 7.5 GHz in these ob¬ 
servations. The spectral resolution was 0.488 MHz per 
channel. J1924-292, J0132-169, Uranus (Neptune for 
band 3 observations) were used for bandpass, phase, and 
flux calibrations. Detailed observational parameters are 
shown in Table 1. 

We used the delivered calibrated data and mapping 
was accomplished using the clean task in CASA (Mc- 
Mullin et al. 2007). We made the data cubes with a 
velocity width of 5 km s“^ for the ^^CO line and 30 km 
s“^ for the other lines. All maps in this paper, except 
for ^^CO (3-2), are reconstructed with a Briggs weight¬ 
ing (robust = 0.5; Briggs & Cornwell 1992) and analyzed 
with MIRIAD and AIPS. The ^^CO (3-2) images are cre¬ 
ated with uniform weighting (see §3.2). The synthesized 
beam size of the (1-0), i®CO (1-0), (3-2), 

and HCN (4-3) were 1"97 x 1"35 (P.A. = 82.3 deg.), 
T.77 X 1"20 (P.A. = 85.8 deg.), 1"64 x 1"17 (P.A. = 
112.6 deg.), and 0746 x 0'.'38 (P.A. = 51.5 deg.), respec¬ 
tively. We also detected CN (l 3 / 2 ~ 0 i/ 2 ), CN (I 1 / 2 -O 1 / 2 ), 
CS (2-1), CH 3 OH (2fc-lfc), and CS (7-6) line emission for 
the first time in VV 114. The properties of these molec¬ 
ular lines are summarized in Table 2. All images which 
we constructed are corrected for primary beam attenu¬ 
ation. The on-source times of band 3 and band 7 were 
about 40 minutes and 80 minutes, and the rms noise 
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Fig. 2.— (a) (1-0) integrated intensity image overlaid on the HST/ACS/F435W image of VV 114. The contours are 0.2, 0.4, 0.8, 

1.6, 3.2, 6.4, 12.8, 25.6, and 33.0 Jy km s“^. The dashed green box shows an imaging field of other lines and continuum of this work 
except for the CO and lines, (b) (1-0) velocity field image. The velocity field image in color scale ranges from 5600 km s~^ 

to 6200 km s“^. The dashed black lines represent tidal arms of VV 114. The dashed red line tracks the filamentary structure detected 
in images of other lines and dust continuum, and the dashed circle shows the overlap region, (c) (1-0) velocity dispersion image. 

The velocity dispersion image in color scale ranges from 0 km s“^ to 120 km s“^. (d) The same as (a) but for (1-0). The contours 

are 0.02, 0.04, 0.08, 0.16, 0.32, and 0.64 Jy km s“^. (e/f) The same as (b/c), respectively, but for (1—0). (g) The same as (a) but 

for CS (2-1). The contours are 0.04, 0.08, 0.16, and 0.28 Jy km s~^. (h/i) The same as (b/c), respectively, but for CS (2-1). The beam 
size of each line is shown in the bottom-left of the images (Table 2). The red crosses show the positions of the nuclei defined by the peak 
positions of the Ks-band observation (Tateuchi et al. 2012). 


levels of the channel maps with 30 km s“^ resolution 
are 1.0 mJy beam“^ and 0.8 mJy beam“^, respectively. 
Furthermore, we made continuum maps at each observ¬ 
ing frequency by adding the line-free channels. The rms 
level of the continuum images were 0.05 mJy beam“^, 
0.11 mJy beam“^, and 0.07 mJy beam“^ for band 3, 
band 7 in the compact conhguration, and band 7 in 
the extended configuration, respectively. The continuum 
emission was subtracted in the rtti-plane before making 
the line images. Throughout this paper, the pixel scales 
of the band 3 and the band 7 images are set to 0'.'3/pixel 
and O'.'08/pixel, respectively, and only the statistical er¬ 
ror is considered unless mentioned otherwise. The sys¬ 
tematic error on the absolute flux is estimated to be 
5% and 10% for both sidebands in band 3 and band 
7, respectively. 

In the following sections, we estimate the missing flux 
of each molecular line for which the single dish data are 


available in literature. Although the effect of missing 
flux becomes critical when we evaluate the global gas 
properties and the corresponding line ratios, the effect 
is negligible when we discuss structures that are smaller 
than the “maximum recoverable scale” (MRS) of each 
conhguration of ALMA. This is estimated from the min¬ 
imum baseline lengths of the assigned antenna conhgu- 
rations and the observed frequencies. The MRS of our 
observations are 8" and ~ 7" in band 3 and band 7, 
respectively (Table 1). Therefore the missing hux effect 
in this paper is negligible, since we derive physical pa¬ 
rameters (e.g., molecular gas mass) only for structures 
smaller than ^ 2". 

3. RESULTS 

Molecular line and continuum images are shown in Fig¬ 
ures 2, 3, 4, 5, and 6. The channel maps and the spectra 
of all line emissions are shown in Appendix A.2 and A.3. 
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Fig. 3.— The same as Figure 2 but for (a, b, and c) CH 3 OH (2^-1^), (d, e, and f) CN (li/ 2 “ 0 i/ 2 )) (§> h, and i) CN (l 3 / 2 “ 0 i/ 2 )- 

(a) The contours are 0.02, 0.04, 0.08, 0.16, and 0.32 Jy km s“^ (d) The contours are 0.04, 0.08, 0.16, 0.32, and 0.50 Jy km s“^. (g) The 
contours are 0.04, 0.08, 0.16, 0.32, 0.64, 1.00, and 1.20 Jy km s“^. 


3 . 1 . Line Emissions in Band 3 
3.1.1. (1~0) 

The integrated intensity, velocity field, and velocity 
dispersion maps of VV 114 are shown in Figures 2a, 
2b, and 2c, respectively. The total (1-0) inte¬ 

grated intensity of VV 114 is 594.6 ± 1.6 Jy km s“^, 
which is 1.3 times larger than that detected using the 
NRAO 12 m telescope (461 Jy km s“^; Sanders et al. 
1991). This is because the pointing center for the NRAO 
12 m observation was 25"0 southwest of the CO cen¬ 
troid identified from the ALMA map (NRAO 12 m; 
01h07m45.7s, -17d30m36.5s; CO centroid; 01h07m47.2s, 
-17d30m25.8s). At the adopted distance of VV 114 
(86 Mpc), the 1'.'97 x 1"35 beam of the (1-0) ob¬ 
servation gives us a resolution of 790 pc x 540 pc. The 
two crosses shown in all images represent the peaks ob¬ 
tained from the miniTAO/ANIR ATs-band observation, 
and we regard them as the progenitor’s nuclei. 

The integrated (1-0) intensity map of VV 114 

(Figure 2a) shows that the diffuse/cold gas forms two 
arm-like structures and a filamentary structure located 
at the center of the image. The global gas distribution 
is consistent with the previous (1-0) observations 


(Yun et al. 1994). The southeastern (SE) arm clearly 
follows the tidal arm seen in the HST/ACS image (Fig¬ 
ure 1; Evans 2008), while the northwestern (NW) arm 
has no counterpart in any other wavelengths. The re¬ 
gion from the center of VV 114 to the eastern nucleus 
shows a strong concentration of molecular gas (~ 5" 5 
west of the eastern nucleus), and we refer to this region 
as the “overlap” region with a molecular “filament” (see 
Figure 2). 

The (1-0) velocity field map of VV 114 (Fig¬ 

ure 2b) shows a significantly broad velocity range across 
the galaxy disks (~ 600 km s“^). The SE arm has a blue- 
shifted velocity from 5650 km s“^ to 5920 km s“^, while 
the NW arm has a red-shifted velocity from 5950 km s“^ 
to 6160 km s“^. One possibility for the larger veloc¬ 
ity width in the SE arm may be a highly inclined tidal 
arm. Two other arm-like features are also detected in 
the (1-0) observations. One arm is located ~ 4"0 

northeast of the eastern nucleus and shows an arc around 
the eastern nucleus in the velocity range of 5810 km s“^ 
to 6180 km s“^. The other arm is located ~ 10"0 
west of the SE arm and has a strong peak (~ 262.5 ± 
1.0 Jy km s“^) in the velocity range of 5610 km to 
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Fig. 4.— The same as Figure 2 but for (a, b, and c) (3-2), (d, e, and f) HCN (4—3), and (g, h, and i) HCO"*" (4-3). (a) The 

contours are 2, 4, 8, 16, 32, 64, 128, and 170 Jy km s“^ (d) The contours are 0.04, 0.08, 0.16, 0.32, 0.64, 1.28, and 1.80 Jy km s“^. (g) 
The contours are 0.04, 0.08, 0.16, 0.32, 0.64, 1.28, and 2.40 Jy km s“^. 


5900 km s-t 

The overlap region has the highest velocity dispersion 
(~ 110 km s“^) (Figure 2c). The NW arm has an average 
velocity dispersion of ~ 30 km s“^, while the SE arm has 
~ 40 km s“^. These values are signihcantly higher than 
the dispersions seen in Giant Molecular Clouds (GMCs) 
in the LMC (2 - 14 km s“^; Minamidani et al. 2008; Fujii 
et al. 2014) and slightly higher than that in Giant Molec¬ 
ular Associations (GMAs) in the Antennae galaxy (6 - 
36 km s“^; Ueda et al. 2012). We suggest that the main 
contribution to the ^^CO ( 1 - 0 ) velocity dispersion is in¬ 
ter cloud turbulent medium along the tidal arm, and/or 
shocked region induced by the tidal interaction, rather 
than the velocity dispersion of the GMCs/GMAs. 

3 . 1 . 2 . ^^CO(1~0) 

The integrated intensity, velocity field, and velocity 
dispersion maps of (1-0) are shown in Figures 2d, 
2e, and 2f, respectively. The integrated (1-0) in¬ 

tensity map of VV 114 (Figure 2d) shows a filamentary 
structure across the galaxy disks, which is consistent 
with the region where the ( 1 - 0 ) filament is de¬ 

tected. The total ^^CO (1-0) integrated intensity is 5.9 


±0.4 Jy km s“^. The strongest peak is located ~ 4"2 
southwest of the eastern nucleus. The ^^CO (1-0) ve¬ 
locity field map of VV 114 (Figure 2e) shows a narrower 
velocity range (5670 - 6000 km s“^) than that of the 
(1-0) emission (5600 - 6200 km s“^). This sug¬ 
gests that the ^^CO (1-0) emission mainly comes from 
two components, the eastern galaxy and the blue-shifted 
component of the overlap region. The ^^CO (1-0) veloc¬ 
ity dispersion map of VV 114 (Figure 2f) shows a high 
velocity dispersion component (^ 100 km s“^) between 
the eastern nucleus and the overlap region. This signif¬ 
icant velocity dispersion may be caused by a superpo¬ 
sition of clouds (see the double-peak spectrum at R39 
shown in Appendix A. 3). 

3.1.3. CS (2-1) and CH 3 OH (2k-lk) 

The GS (2-1) and CH3OH (2fc-lfc) lines are only de¬ 
tected at the overlap region (Figures 2g, 2h, 2i, 3a, 3b, 
and 3c). This is the first detection of the GH3OH (2fc- 
Ifc) emission in a merger-induced overlap region. We 
observed the blended set of 2i - li (rarest = 96.756 GHz, 
A.p/fc = 28.0 K), 2o - lo A (i.,est = 96.745 GHz, A.p/fc 
= 20.1 K), 2o - lo A+ (zzrest = 96.741 GHz, E^^/k = 
7.0 K), and 2_i - l_i E (rarest = 96.739 GHz, E^p/k = 
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Fig. 5.— (a) 340 GHz continuum flux image of VV 114E. The flux 
in color scale ranges from 0 mjy beam“^ to 3.4 mjy beam“^. (b) 
CS (7-6) integrated intensity image of VV 114E. The flux in color 
scale ranges from 0 Jy km s~^ to 0.45 Jy km s~^. (c) HCN (4- 
3) integrated intensity image of VV 114E. The flux in color scale 
ranges from 0 Jy km s“^ to 2.0 Jy km s“^. (d) HCO+ (4-3) 

integrated intensity image of VV 114E. The flux in color scale 
ranges from 0 Jy km s“^ to 2.8 Jy km s~^. The beam size of 
each line and continuum is shown in the bottom-left of the images 
(Table 2). The red cross shows the position of the eastern nucleus 
defined by the peak position of the Ks-band observation (Tateuchi 
et al. 2012). 

12.5 K), thermal transitions of CH 3 OH (hereafter desig¬ 
nated the 2fc - Ife transition). The distribution of these 
molecular lines is clearly different from the other dense 
gas tracers detected in the current program. The peaks of 
CS (2-1) and CH 3 OH (2^-1^) are coincident with one of 
the peaks of (1-0) to within O'.'S. The total CS (2- 
1) and CH 3 OH (2fe-lfc) integrated intensities are 0.4 ± 
0.1 Jy km s“^ and 0.5 ± 0.1 Jy km s“^, respectively. The 
signal to noise is too low to resolve the velocity structure. 

3.1.4. CN (l3/2~0i/2) cwd CN (li/2~0i/2) 

Two radical CN rotational transitions TV = 1 - 0 (J = 
3/2 - 1/2 and 1/2 - 1/2) are detected at the eastern nu¬ 
cleus. The J = 3/2 - 1/2 transition is extended toward 
the overlap region (Figures 3d, 3e, 3f, 3g, 3h, and 3i). 
We can not resolve their multiplet because of the coarse 
frequency resolution (11.5 MHz ~ 30 km s“^). Because 
the critical density of CN is high (^ 10® cm“^), the CN 
emission mainly comes from denser gas regions than re¬ 
gions traced by ^^CO (1-0). The J = 3/2 - 1/2 transition 
shows a similar distribution to the ^^CO ( 1 - 0 ) emission, 
but it is less extended over the overlap region. The total 
CN (li/ 2 " 0 i/ 2 ) and CN (l 3 / 2 " 0 i/ 2 ) integrated intensities 
are 2.0 ± 0.1 Jy km s“^ and 5.4 ± 0.3 Jy km s“^, respec¬ 
tively. The highest velocity dispersion in the CN (I 3 / 2 - 
O 1 / 2 ) image is also detected between the eastern nucleus 
and the overlap region, and this is likely caused by a su¬ 
perposition of clouds similar to the case of the ^^CO ( 1 - 0 ) 
image (see Appendix A. 3). 

3.2. Line Emission in Band 7 


3.2.1. CO (3-2) 

The ^^CO (3-2) emission maps are presented in Fig¬ 
ure 4. The estimated missing flux in our ALMA ob¬ 
servation is 21 ± 1 % {James Clerk Maxwell Telescope 
(JCMT): 2956 ± 133 Jy km s-^ and ALMA: 2343.7 ± 
4.7 Jy km s“^; Wilson et al. 2008; Saito et al. 2013). Al¬ 
though our ^^CO (3-2) observation recovers more flux 
than the Submillimeter Array (SMA) observation (1530 
± 16 Jy km s“^; the missing flux = 48 ± 15 %; Wil¬ 
son et al. 2008), there are significant negative sidelobes 
at the north and south of the image which is likely the 
cause of missing flux. We made the CLEANed image 
with a uniform uv weighting to minimized the sidelobe 
level (Thompson et al. 2001). 

The ^^CO (3-2) integrated intensity map of VV 114 
(Figure 4a) shows two arm-like structures and a filamen¬ 
tary structure similar to the ^^CO ( 1 - 0 ) image, and the 
strongest peak is at ~ 5'.'5 west of the eastern nucleus. 
The global gas distribution is consistent with the previ¬ 
ous ^^CO (3-2) observations (lono et al. 2004; Wilson et 
al. 2008). The ^^CO (3-2) velocity held map of VV 114 
(Figure 4b) also shows signihcant broad velocity range 
across the galaxy disks (~ 600 km s“^), similar to the 
^^CO (1-0) velocity held map. The SE arm has a blue- 
shifted velocity from 5650 km s“^ to 5920 km s“^, while 
the NW arm has a red-shifted velocity from 5950 km s“^ 
to 6160 km s“^. Other two arm-like features are also 
detected. One located ~ 4'.'0 northeast of the eastern 
nucleus shows an arc around the eastern nucleus and has 
red-shifted velocities from 5810 km s“^ to 6180 km s“^. 
This arm coincides with the NE arm detected in the 
^^CO (1-0). The other one located at ~ lO'.'O west of the 
SE arm has a strong peak (~ 262.5 ± 0.9 Jy km s“^) and 
blue-shifted velocities from 5610 km s“^ to 5900 km s“^. 
This arm also coincide with the SW arm detected in the 
(1-0). From the (3-2) velocity dispersion 

map of VV 114 (Figure 4c), we find that the overlap re¬ 
gion has the highest velocity dispersion (~ 110 km s“^). 
The velocity dispersion of the NW arm is ~ 30 km s“^, 
while the SE arm is ~ 60 km s“^. 

3.2.2. HCN (4-3) and HCO+ (4-3) 

The HCN (4-3) and HCO’*' (4-3) images are shown in 
Eigures 4d, 4e, 4f, 4g, 4h, and 4i. While the HCN (4-3) 
emission is only seen near the eastern nucleus of VV 114 
and is resolved into four peaks, the HCO’*' (4-3) emis¬ 
sion is more extended and has at least 10 peaks in the 
integrated intensity map. The total integrated intensi¬ 
ties of HCO+ (4-3) and HCN (4-3) are 15.3 ± 0.4 Jy km 
s“^ and 4.4 ± 0.2 Jy km s“^, respectively. The higher 
HCO’*' (4-3) flux observed with the SMA (17 ± 2 mJy, 
Wilson et al. 2008) using a 278 x 270 beam is likely at¬ 
tributed to missing flux by the ALMA observation. A 
compact component in the eastern nucleus is unresolved 
with the current resolution, and the upper limit to the 
size is 200 pc. The HCN (4-3) emission is not detected in 
the overlap region, where both the high ^^CO ( 1 - 0 ) ve¬ 
locity dispersion and the significant CH 3 OH (2fc-lfc) and 
HCO+ (4-3) detection suggest the presence of shocked 
gas (Krips et al. 2008). We concluded in paper I from 
their source size, line widths, and the relative strengths of 
HCN (4-3) and HCO^ (4-3) that the unresolved eastern 
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Fig. 6 .— (a) The 110 GHz continuum flux image overlaid on the HST/ACS/F435W image of VV 114. The contours are 0.10, 0.20, 0.40, 
0.80, and 1.60 mjy bearn”^. The red crosses show the positions of the nuclei defined by the peak positions of the Ks-band observation 
(Tateuchi et ah 2012). (b) The low resolution 340 GHz continuum flux image overlaid on the HST/ACS/F435W image of VV 114. The 
contours are 0.22, 0.44, 0.88, 1.76, 3.52, and 7.04 mJy bearn”^. A strong point source at the eastern edge of the image is a distant 
star-forming galaxy, ALMA J010748.3 - 173028 (see Tamura et ah 2014). (c) The high resolution 340 GHz continuum flux image overlaid 
on the HST/ACS/F435W image of VV 114. The contours are 0.14, 0.28, 0.56, 1.12, and 2.24 mJy beam“^. The red cross shows the 
position of the eastern nucleus defined by the peak position of the Ks-band observation (Tateuchi et ah 2012). The beam size of each 
continuum is shown in the bottom-left of the images (Table 2). 


nucleus harbors an obscured AGN, and the dense clumps 
in the western galaxy are related to extended starbursts. 

3.2.3. CS (7-6) 

The CS (7-6) emission has the highest critical density 
(ricr — 10^ cm“^) of all of the lines detected in our obser¬ 
vations. The CS (7-6) emission is marginally (S/N ~ 4) 
detected at the eastern nucleus (Figure 5), and the total 
flux is 0.5 ± 0.1 Jy km s“^. 

3.3. Continuum Emission 

The continuum image at 110 GHz shows a filamentary 
structure similar to the molecular line image (Figure 6a). 
We construct low resolution (1733 x 1712) and high res¬ 
olution (0745 X 0738) images of the 340 GHz continuum 
(Figures 6b, and 6c) using the combined data (compact -|- 
extended) and the extended configuration data, respec¬ 
tively. We find that the filamentary structure and the 
unresolved eastern nucleus are both present in dust con¬ 
tinuum. The total flux of the 110 GHz and the low reso¬ 
lution 340 GHz continuum emission are 10.3 ± 0.2 mJy 
and 38.6 ± 0.3 mJy, respectively. The estimated missing 
flux relative to the JCMT 340 GHz observation (Wilson 
et al. 2008) is 75 ± 4 % (SMA: 79 ± 7 %). The difference 
in the recovered flux between (3-2) and 340 GHz 

continuum emission may be caused by the difference in 
the distribution. The 110 GHz and 340 GHz continuum 
emission is detected at the eastern nucleus (S/N ^ 50 
and 70) and the filamentary structure (S/N ^ 8 and 24) 
identified in the ^^GO (1-0) image, both with high sig¬ 
nificance. 

4. SPATIALLY RESOLVED LINE RATIOS 

We assign 39 “R” boxes (270 x 270; R1 - R39; see 
Figure 7) for the band 3 and (3-2) data and 15 

smaller “S” boxes (1"2 x 172; SO - S14; see Figure 8) for 
the rest of the data to estimate the physical parameters, 
such as the molecular gas mass (Mhj), dense gas mass 
(Mdense), dust mass (Mdust), star formation rate (SFR), 
kinetic temperature (Tkin), gas density (nHj), gas col¬ 
umn density (V(H 2 )), and molecular abundance relative 


to H 2 ([A]/[H 2 ]). The positions of the boxes are cho¬ 
sen to cover the CO (3-2) emission (R1 - R39) and the 
HCO+ (4-3) emission (SO - S14). The sizes of the boxes 
are chosen such that they are comparable to the beam 
size. Before deriving the parameters and line ratios at 
each box, we first matched the uv range between our 
data set and reconstructed the integrated intensity im¬ 
age of each line. The shortest baseline lengths are set to 
13.5 kA and 40.0 kA for the molecular lines in the band 
3 and the band 7, respectively, and the images are con¬ 
volved into the same resolution (270 x 175 with a P.A. 
of 83 deg, 172 x T.'O with a P.A. of 119 deg). For each 
ratio, the two integrated intensity images were expressed 
in the units of K km s“^ before calculating the ratio at 
locations where both lines are detected above 3 a. The 
derived box-summed spectra are listed in Appendix A.3. 
We carried out a multi Gaussian fit (one - three compo¬ 
nents) to reproduce the box-summed spectra, and labeled 
the components as “a”, “b”, and “c” from the bluest peak 
(e.g., the bluest peak at R21 is labeled as R21a). 

4.1. (3-2)/^^CO (1-0), R3_2Ii-o 

The (3-2)/^^CO (1-0) ratio, i? 3 _ 2 /i_o, can be 

used as an indicator of the dense/warm gas content rela¬ 
tive to the total molecular gas. The R 3 - 2 / 1-0 of VV 114 
varies from 0.2 to 0.8, as shown in Figure 7 (left) and 
Table 6. This range is larger than the same ratios de¬ 
rived for normal spirals, which is typically 0.15 - 0.5 
when observed with a similar linear resolution (Warren 
et al. 2010). At the edge and the center of the filament, 
i? 3 _ 2 /i-o is higher (0.53 - 0.69) than the highest peaks 
of each arm (^ 0.4). This suggests that the GO emit¬ 
ting gas at the filament have higher excitation conditions 
than normal spirals, while the conditions of each arm of 
VV 114 are consistent with arms and nuclei of normal 
spirals. The i? 3 _ 2 /i-o at the eastern nucleus is 0.76 ± 
0.01. It is suggested that the i? 3 _ 2 /i_o is much higher 
(3.12 ± 0.03 in NGG 1068; Tsai et al. 2012) for gas sur¬ 
rounding an AGN, and the low i? 3 _ 2 /i_o in VV 114 may 
be due to the difference in filling factor (160 x 140 pc 
beam averaging for NGG 1068, while 800 pc box averag¬ 
ing for VV 114 ). It is possible, however, that the nuclear 
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Fig. 7.— (a) The Rs— 2 / 1—0 image. The ratio in color scale ranges from 0 to 1. The white crosses show the positions of the nuclei 
defined by the peak positions of the Ks-band observation (Tateuchi et al. 2012). (b) The R 12/13 image. The ratio in color scale ranges 
from 0 to 40. (c) Locations of 39 boxes (R1 — R39) that are used to calculate the line ratios and physical parameters. For each ratio, the 
two integrated intensity images were convolved to the same resolution and expressed in units of K km s~^ before calculating the ratio at 
locations where both lines are detected above 3 cr. The black crosses show the positions of the nuclei defined by the peak positions of the 
Ks-band observation (Tateuchi et al. 2012). The beam size of each line ratio is shown in the bottom-left of the images. 


excitation conditions are different from source to source. 

4 . 2 . (l-O)A^CO (1-0), Ri2/i3 

In general, the lines has higher optical depths 

than the (1-0) line. Therefore, the measured 

(1-0)/^^CO (1-0) line intensity ratio, i?i 2 /i 3 , gives 
a lower limit to the CO/^^CO abundance ratio (hereafter 
[CO]/[^^CO]). We present the R 12/13 image of VV 114 
in Figure 7 (center). The i?i 2 /i 3 increases from the 
arms (< 17) to the filament (15 - 32). Observationally, 
7 ?i 2/13 increases towards the central region of galaxies 
(Aalto et al. 1995), where the gas is generally warmer 
and denser. Aalto et al. (1995) suggest that the moder¬ 
ate optical depth of ( 1 - 0 ) emission and/or the high 
[CO]/[^^CO] environment can increase the i?i 2 /i 3 in nu¬ 
clei of U/LIRGs. In order to understand which of the 
two (optical depths or abundances) is dominant, we cal¬ 
culated and mapped the optical depth of the ( 1 - 0 ) 
and the (1-0) as shown in Table 14 and Figure 10. 
We provide an interpretation of these results in §5.2. 

4.3. HCN (4-3)/HCO+ (4-3), i?HCN/HCO+ 

In paper I, the HCN (4-3) and HCO''" (4-3) maps of 
VV 114 allowed us to investigate the central region at 
200 pc resolution for the first time, and we find that 
both the HCN (4-3) and HCO’*' (4-3) in the eastern nu¬ 
cleus are compact (< 200 pc), and broad [290 km s“^ 
for HCN (4-3)]. We present the HCN (4-3)/HCO+ (4- 
3), Ahcn/hco+: image of VV 114 in Figure 8 . From 
the higher i?HCN/HCO+ along with the past X-ray and 
NIR observations, we suggest the presence of an obscured 
AGN in the eastern nucleus. We also detect a 3 - 4 kpc 
long filament of dense gas, which is likely to be trac¬ 
ing the active star formation triggered by the ongoing 
merger, and this is consistent with the results from the 
numerical model by Teyssier et al. (2010) who predict 
that the fragmentation and turbulent motion of dense 
gas across the merging disk is responsible for forming 
dense gas clumps with masses of 10 ® - 10 ® M 0 . 

We present the i?HCN/HCO+ image in Figure 8. The 
overlap region does not show significant HCN (4-3) emis¬ 


sion, and we provide the 3 a upper limit in Table 7. Three 
out of the four boxes (i.e., SI - S3) in the eastern nu¬ 
cleus have low i?HCN/HCO+ (< 0.5) whereas SO has a high 
i?HCN/HCO+ (1-34 ± 0.09). It is suggested that such a 
high value is only produced around AGN environments 
(e.g., Kohno et al. 2001; Harada et al. 2013; lono et al. 
2013; Izumi et al. 2013; Imanishi & Nakanishi 2014). 

5. DERIVATION OF PHYSICAL PARAMETERS 

In this section, we derive the molecular gas mass (§5.1), 
and the physical parameters using the radiative transfer 
code RADEX (§5.2) for each box defined in §4. The 
column density is derived using the optically thin ^®CO 
line under the LTE assumption. We estimate the beam 
filling factor <I>a and the relative molecular abundance 
of molecule X (hereafter expressed as [A']/[H 2 ]) (§5.3). 
Finally, we calculate the dust mass using the 340 GHz 
continuum emission (§5.4). 

5.1. Molecular Gas Mass Derivation 

The molecular gas mass Mx is derived by; 

A/x = ax Ax [Mq], (1) 

where ax is the molecular line luminosity-to-H 2 mass 
conversion factor and A'x is the velocity integrated 
flux (Solomon & Vanden Bout 2005). We use 
the conversion factor known to be appropriate for 
U/LIRGs (aco = 0.8 Mq (K km s“^ pc^)“^; Downes 
& Solomon 1998). This is consistent with the value 
derived by Sliwa et al. (2013) in VV 114 (aco = 
O.hlg'g M 0 (K km s“^ pc^)“^). The molecular gas mass 
derived at the boxes defined in §4 ranges between 0.2 
X 10® (ff) and 4.8 x 10® (ff) M 0 (Table 10). We 
also calculate the dense gas mass M^ense using ancN = 
10 M 0 (K km s“^ pc^)“^ (Gao & Solomon 2004) and the 
HCN (4-3) luminosity which is converted to the HCN (1- 
0) luminosity using HCN (4-3)/HCN (1-0) = 0.63 (pa¬ 
per I; Imanishi et al. 2007). The dense gas mass ranges 
between 1.8 x 10® and 3.8 x 10^ M 0 

(Table 11). 

We note that the CO luminosity-to-H 2 mass conversion 
factor, aco I is very uncertain, and varies signihcantly 










Molecular gas in VV 114 


9 




RA (JSOOO) 



jV4'f‘6 irs Af.\ 47'0 4«*9 


RA (J200a} 


Fig. 8. — (a) The ^hcn/hco+ The ratio in color scale ranges from 0 to 1. The white cross shows the position of the eastern 

nucleus defined by the peak position of the Ks-band observation (Tateuchi et ah 2012). (b) The ^hcn/hco+ near the nucleus of 

VV 114E. The ratio in color scale also ranges from 0 to 1. (c) Locations of 15 boxes (SO - S14) that are used to calculate the line ratios 
and physical parameters. For each ratio, the two integrated intensity images were convolved to the same resolution and expressed in units 
of K km s~^ before calculating the ratio at locations where both lines are detected above 3 a. The black cross shows the position of the 
eastern nucleus defined by the peak position of the Ks-band observation (Tateuchi et al. 2012). The green open squares are the “R” boxes 
shown in Figure 7c. 


from source to source (0.4 - 0.8 for LlRGs; Downes & 
Solomon 1998; Yao et al. 2003; Papadopoulos et al. 2012; 
Bolatto et al. 2013). It may be possible that aco varies 
from region to region within a galaxy. While one would 
ideally adopt a spatially varying aco for a better quan¬ 
tification of the H 2 mass, such a study is beyond the 
scope of this present paper. For simplicity, here we adopt 
a constant aco across all regions in VV 114, bearing in 
mind that the uncertainties could be as large as a fac¬ 
tor of two. The same applies to ancN (Gao & Solomon 
2004). 

5.2. Radiative Transfer Analysis using RADEX 

We used the non-LTE radiative transfer code RADEX 
(van der Tak et al. 2007) and varied the parameters un¬ 
til the residuals between the observed line fluxes and the 
modeled line fluxes are minimized in a sense. We as¬ 
sumed a uniform spherical geometry [dv = 1.0 km s“^), 
and derived the physical conditions of molecular gas 
(Tkin, and 7V(H2)). RADEX uses an escape prob¬ 
ability approximation to solve the non-LTE excitation 
assuming that all lines are from the same region. Since 
the molecular lines in the band 7 have significantly higher 
critical densities than that in the band 3, we used two sets 
of molecular lines; (case 1) 2"0 box-summed ^^GO (1- 
0), ^^GO (1-0), and ^^GO (3-2), and (case 2) 1"2 box- 
summed HGN (4-3), HCO+ (4-3), ^^GO (3-2), and 
^^GO (1-0), to solve for the degeneracy of the physi¬ 
cal parameters. In case 2, we made the uv and beam- 
matched HGN (4-3), HGO+ (4-3), and (3-2) im¬ 
ages (1"2 X I'.'O resolution with the P.A. = 119 deg.), 
and we defined three HCO+ (4-3) peaks as EO, El, and 
E2 (Figure 9). We also use the uv and beam-matched 
^^GO (1-0) data to constrain the 7V(H2), allowing us to 
vary the [HGN]/[HGO'''] in case 2. All line parameters, 
such as the upper state energies and the Einstein coef- 
hcients, were taken from the Leiden Atomic and Molec¬ 
ular Database (LAMDA; Schoier et al. 2005). In order 
to find the set of physical parameters that can repro¬ 
duce the observed line intensities, we run RADEX by 
varying Tkin, ?^H 2 j and 1 V(H 2 ) for case 1, and Tkm, ?^H 2 i 
and [HGN]/[HCO+] for case 2. The adopted iV(H 2 ) are 
10^^’^, 10^^ ®, and 10^^ ® cm“^, at EO, El, and E2, re¬ 


spectively. 

We varied the gas kinetic temperature within a range 
of Tkin = 5 - 300 K using steps of dTkin = 5 K, and a 
gas density of nH 2 = 10^ - 10^ cm“^ using steps of dnH 2 
= 10°'^ cm“^. For case 1, we fixed [^^GO]/[H 2 ] = 1.4 x 
10“® (Davis et al. 2013) and [GO]/[^^GO] = 70, which 
are the Galactic values (Wilson & Rood 1994). In case 
2, we changed the parameters, Tkin = 5 - 400 K using 
steps of dTkin = 5 K, nH 2 = 10^ - 10^ cm“^ using steps 
of dnH 2 = 10°'^ cm“^, and fixed [CO]/[H 2 ] = 1.0 x 10“^ 
and [HC 0 +]/[H 2 ] = 1.0 x 10“®, which are the standard 
values observed in Galactic molecular clouds (Blake et al. 
1987). We varied [HCN]/[HCO'''] from 1 - 10, in steps of 
one. The parameters we used are summarized in Table 9. 
We list the results that are within the 95 % confidence 
level with 3-degree of freedom (x^ < 7.81) (Tables 12 
and 13). Finally, we created velocity-averaged channel 
maps of nH 2 and the optical depth of the transitions (Fig¬ 
ure 10). 

We note that the uncertainty of the V(H 2 ) for case 
2 did not strongly affect the results, while that of the 
[GO]/[^^GO] for case 1 changed. The effect of varying the 
[GO]/[^^GO] will be discussed in §5.2.1. Future multi¬ 
transition HCN/HCO+/GO/13GO imaging will help us 
to derive these parameters directly. 

5.2.1. Case 1 

The (box-averaged) kinetic temperature near the east¬ 
ern nucleus (R21a) is constrained to within 25 - 90 K (the 
best ht is 50 K), as shown in Table 12. The Tex (58.8 
± 2.9 K) obtained from the LTE assumption at R21a 
(see §5.3) is higher than the best-fitted Tkin- In fact, 
we also find five regions (RlOb, Rllb, R14, R16, and 
R25a) that show similarly high excitation temperatures. 
Four out of five regions are in the central filament. In 
general, spontaneous emission dominates over collisional 
excitation in sub-thermally excited conditions, and hence 
Tex should be lower than Tkin- One reason for this dis¬ 
crepancy could be attributed to the incorrect assump¬ 
tion of [GO]/[^^GO]. By varying this abundance ratio, 
we find that the temperature reversal (i.e. Tkin > Tgx) 
occurs only when [CO]/[^^CO] > 150. This is consistent 
with the results obtained by Sliwa et al. (2013) who used 
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Fig. 9.— (top) The uv- and beam-matched (grey color) (3-2), (red contour) HCN (4-3), and (blue contour) HCO"*" (4-3) images. 

The integrated intensity of the (3-2) in color scale ranges from 0 Jy km s~^ to 100 Jy km s“^. The contours are 5, 10, 20, 40, 

and 50 Jy km s“^ for HCN (4—3), and 12, 24, 48, and 96 Jy km s“^ for HCO"*" (4—3). The white cross shows the position of the eastern 

nucleus defined by the peak position of the Ks-band observation (Tateuchi et al. 2012). (bottom) 1'.'2 box-summed spectra of (black line) 

(3-2), (red dashed line) HCN (4-3) x 10, and (blue dashed line) HCO"t (4-3) x 10 at the each box, labeled EO — E2. The spectra 

are taken from the ALMA data cubes after correcting the cubes for the primary beam attenuation and convolving them to 1^2 x F.'O 
resolution (P.A. = 119 deg.). 


RADEX along with their multi CO and line data 

to find evidence of a cold/dense molecular gas compo¬ 
nent with extremely high [CO]/[^^CO] of 229, which is 3 
times higher than that of the Galactic value (Wilson & 
Rood 1994). 

The derived Tkm at the other regions are generally 
higher than 100 K. The derived Tex in each region are 
typically 10-40 K, which may suggest sub-thermal con¬ 
ditions. The kinetic temperatures derived at the SE and 
NW arms are estimated to be < 90 K, with higher tem¬ 
perature at the NW arm. The NW arm is also associated 
with relatively strong Paa emission and ATs-band emis¬ 
sion, which is consistent with the higher relative temper¬ 
ature due to star-forming activities (Minamidani et al. 
2008). However, this is inconsistent with the general un¬ 
derstanding that strong tidal shear in tidal arms prevents 
active star formation to occur (Aalto et al. 2010). 

The derived nua in most of the boxes are less than 
10 ^ ° cm“^, which is consistent with the critical densi¬ 
ties of the low-J CO lines observed here. The highest 
density of 10^ "^ - 10® ° cm“® is estimated at R21a, and 
this is consistent with the location of the eastern nu¬ 
cleus. Since we also observed the strongest HCN (4-3) 
and HCO’*' (4-3) emission at R21a at the same line-of- 
sight velocity (lono et al. 2013, see also Appendix A. 3), 
it is possible that the main contribution to the CO emis¬ 
sion at R21a arises from dense gas (10®'^ - 10® ° cm“°) 
near the eastern nucleus, with a minor contribution from 
the diffuse gas clouds along the same line of sight ob¬ 
served within the same beam. In contrast to the eastern 


nucleus, the boxes that cover the western galaxy (R1 - 
Rll and R26 - R29) show moderately dense condition of 
102 ° _ o cm“°. This extended and moderately dense 
gas is associated with the disk-like structure seen in op¬ 
tical images (Evans 2008), and the star formation traced 
in Paa emission and UV/X-ray emission (Grimes et al. 
2006; Tateuchi et al. 2012). We note that the strongest 
off-nuclear Paa peak (R2'7 in Table 10; SFR = 3.15 ± 
0.05 Mq yr“^) coincides with relatively low gas density 
(~ 10® ° cm“®). The density of the surrounding region 
labeled R25a is similar (10® ® - 10® ° cm“®) and this is 
comparable to the nucleus of the eastern galaxy. The 
secondary Paa peak (R29; SFR = 0.92 ± 0.05 Mq yr“®) 
is not associated with any molecular line emission. 

It is usually believed that the ®^CO ( 1 - 0 ) emission is 
optically thick (rco 1), while the ®®CO (1-0) emission 
is optically thin {tisqq ^ 1 ) even in luminous mergers 
(Davis et al. 2013). In most regions, we find that the op¬ 
tical depth of the ^^CO (1-0) line is ^ 1 (Figure 10). In 
contrast, the ®^CO ( 1 - 0 ) opacity at the eastern nucleus 
and the filament is moderately optically thick (rco !)■ 
However, the elevated T 12/13 at the eastern nucleus (see 
§4) cannot be explained by the relatively low ®^CO (1-0) 
opacity alone (the opacity has to be rco 0 . 1 ; see also 
Wilson et al. (2009)). Finally, we find that indeed the 
^®CO ( 1 - 0 ) emission is optically thin (risco 1 ) aver¬ 
aged over the whole galaxy, except for the southern dust 
lane (risco = 0-3 - 1-5). 

From these results, we suggest that the peak of the 
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Fig. 10.— The channel maps of the box-averaged RADEX modeling based on the (1-0), (1-0), and (3-2). (top) The 

best fitted values of logarithmic gas density of the (1-0) emission. The value in color scale ranges from 2.0 to 4.0 cm“®. (bottom) 

The best fitted values of optical depth of the (1-0) emission. The value in color scale ranges from 0 to 9. The black crosses show 

the positions of the nuclei defined by the peak positions of the Ks-band observation (Tateuchi et al. 2012). The open squares are regions 
which we cannot solve the RADEX calculations because of non-detection of the (1—0), (1-0), or (3-2) emission. 


molecular gas in the central 800 pc of the eastern galaxy 
is cold (Tkin = 25 - 90 K), dense (riHa = 10^ "^ - 
10 ® ° cm“®), and moderately optically thick (tco(i-o) 
^ 3), while peaks in the overlap region are warm (Tkin 
> 50 K, best-fitted Tkin is 95 and 175 K at R39a and 
R39b, respectively), moderately dense (uhs = 10^ ® - 
10 ^-^ cm“®), and moderately optically thick (tqo(i-o) 

1). The derived density of the eastern galaxy is slightly 
higher than the range of values found in U/LlRGs using 
low-J CO emission with^ kpc resolution (riHa = 10^ ® - 
10^ ® cm“®; Downes & Solomon 1998). In addition, the 
low opacities predicted from these analyses are consistent 
with earlier results that investigate the opacities in M82 
(r = 0.5 - 4.5; Mao et al. 2000) and U/LIRGs (r = 3 - 
10; Downes & Solomon 1998), and the central region of 
NGC 6240 (r = 0.2 - 2; lono et al. 2007). However, the 
derived temperature of the eastern galaxy is inconsistent 
with the high values found in nearby starburst galaxies 
M82, NGC 253, and NGC 6240 (Wild et al. 1992; Jack- 
son et al. 1995; Seaquist & Frayer 2000; lono et al. 2007). 
The disagreement is possibly due to the uncertainties in 
the [CO]/[^®CO], or the difference in the observed molec¬ 
ular gas tracers. 


5.2.2. Case 2 

The values for Tkin, s-iid the optical depth of 

HCN (4-3) and HCO+ (4-3) are shown in Table 8. 
The derived parameters for the unresolved component 
EO, are T^n > 100 K, uhs = 10 ® ° - 10 ®'^ cm“®, and 
[HCNj/pCO+J > 5. The lower limit to the kinetic tem¬ 
perature is higher than those of El and E2, mainly due 
to the unusually high Thcn/hco+ a-^d Rhcn/co- In 


contrast to EO, the derived parameters near El show 
high H 2 densities (riHa = 10^’^ ^ 10^'® cm“®). The over¬ 
lap region (E2), where the star-formation rate (1.70 ± 
0.05 Mq yr“^) is lower than the eastern nucleus, has den¬ 
sities in the range of uhs = 10® ° - 10® ° cm“®. Einally, 
the optical depths for the HCO+ (4-3) and HCN (4-3) 
lines are calculated for each gas clump, yielding thcn — 
0.7 and T{jco+ — 0.2 for EO, thcn — 0.2 and thco+ — 
0.6 for El, and thcn — 0.4 and thco+ — 0-4 for E2. 

The higher linear resolution observations of HCN (4- 
3) and HCO'*' (4-3) toward NGC 1097 (Izumi et al. 
2013) revealed that the gas in the central region of 
NGC 1097 has Tki„ = 70 - 550 K and = 10^-^ - 

10° ° cm“®. Moreover, by comparing to LVG models, 
Krips et al. (2008) found that HCN and HCO"*' emission 
in AGN-dominated sources appears to emerge from re¬ 
gions with lower H 2 densities, higher temperatures, and 
higher HCN abundance relative to starburst-dominated 
(SB-dominated) galaxies. Our results obtained toward 
VV 114 are consistent with these previous results. 

5.3. Filling factor and Column Density under LTE 

In order to determine the bulk properties of the CO 
emitting gas, we used an excitation temperature analysis 
(Davis et al. 2013). The excitation temperature at each 
box can be calculated from 


Tex = Tn In 


($aTo( 1 - e—^co) ' e^o/Tbg _ i 


fb,CO(l-0) 


-1 


-f 1 


(2) 
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where Tq = hv/k [= 5.53 K for (1-0) emission], v 

is the frequency of the transition, h is the Planck’s con¬ 
stant, k is the Boltzmann’s constant, 7b,co(i-o) is the 
brightness temperature of (1-0) emission in Kelvin, 
Tco is the optical depth of the ( 1 - 0 ) emission, 

and Tbg is the cosmic microwave background tempera¬ 
ture (2.73 K). Using Tkin estimated from the RADEX 
calculation (§5.2), we estimate the beam filling factor 
$A, 


_ ^b,CO(l-0) 

Tkin 


(3) 


The optical depth of the (1-0) emission is also 

estimated from the RADEX calculation in §5.2. Assum¬ 
ing that the and CO arise from the same molecular 

cloud, and that the ( 1 - 0 ) is optically thick, we es¬ 

timate the optical depth of a given molecule using, 


since the long wavelength Rayleigh-Jeans (RJ) tail of 
thermal dust emission is often optically thin. In or¬ 
der to compare the Mism with the Mhj (see $5.1) us¬ 
ing spatially-resolved data, we calculated the total ISM 
mass from the 340 GHz continuum emission (Scoville et 
al. 2014). For z^rest ^ 1199 GHz, 




^obs 


,353 GHz') 


3.8 




( 8 ) 


where is the observed flux, Mism is the ISM mass, 
z/obs is the observed frequency, and Erj and Eq are given 

by 


rRj(7d, Z/obs, 2 ) 


hi^ohsil + z)/kTa 

e?ii^obs(i+z)/fc7’d _ I 


(9) 


rx ~ In 



Th,x 

Tb.CO(l-O) 


- 1 ' 


(4) 


where rx is the optical depth of a given transition, and 
2 b,x is the observed brightness temperature for transi¬ 
tion X. Using Tex and rx, we estimate the column density 
for a given molecule from, 


Nx = 


3fc 

+ 1 ) 


«p(^) 

('-“p(-tfe)) 


X 


TX _1_ 

1 - e-x J(Tex) - J(Tbg) 


T^dV 


J{T) 


hv 1 

k exp{hi'/kT) — 1 


(5) 

( 6 ) 


ro = rRj(rd,353 GHz,0). (lo) 

The derived box-summed ISM masses of VV 114 range 
between 5.2 x 10^ and 7.2 x 10® Mq (Table 10). This is 
comparable to the box-summed H 2 masses (Mrj = ( 0.2 
- 4.7) X 10® {^) Mq). We find that the Mism/Mr^ 
ratio is close to unity (0.5 - 2.0, the average Mism/A 7 h 2 
= 0.9 ± 0.1), while the total Mism/Mrs ratio is 0.6 ± 
0.1. This means that the spatially-resolved Mism is a 
good tracer of the “resolved” H 2 mass. However, the 
total Mism underestimates the H 2 mass (even using the 
aco for ULIRGs to derive the Mr^) because the global 
distribution of the 340 GHz continuum emission is sig¬ 
nificantly different from that of the GO (1-0) emission 
(Figures 2 and 6 ). This difference between the 340 GHz 
continuum and the GO (1-0) is also seen in recent obser¬ 
vations of nearby LlRGs (e.g., Sakamoto et al. 2014). 


where ^ is the dipole moment, B is the rotational con¬ 
stant, J is the lower energy level, and f T^dV is the 
integrated intensity (Wilson et al. 2009). The derived 
column densities are listed in Tables 14 and 15. 

5.4. Bust Mass and ISM Mass Derivation from 
340 GHz continuum 

We calculated the dust mass from the 340 GHz 
(880 ^m) continuum emission (Table 3) using (Wilson 
et al. 2008), 

17 

Mdust = 74220 ^340 Dl Mq (7) 

K340 

where £'340 is the 340 GHz flux in Jy and is the lu¬ 
minosity distance in Mpc. We assumed a dust emis- 
sivity, K 340 = 0.9 cm^ g“^, and the dust temperature 
Td of 39.4 K (Wilson et al. 2008). The box-summed 

dust masses ranges between 2.0 x 10 ^ 2.8 x 

10® (Table 10). We note that we used the 

Draine & Lee (1984) dust model for K 340 , because the 
zt 34 o derived from observations has a large error (Hen¬ 
ning et al. 1995). 

Scoville et al. (2014) suggested that the submillimeter 
continuum emission traces the total ISM mass (Mism), 


6 . DISCUSSION 

6 . 1 . Conditions of “Dense” Gas near the Eastern 
Nucleus 

Our RADEX modeling yields lower molecular gas den¬ 
sity near the AGN (rr^ = 10® ° - 10®'^ cm“®) compared 
to the surrounding clumps (10® ® - 10® ° cm“®). Simi¬ 
larly high values are obtained near AGNs in other galax¬ 
ies ( Alonso-Herrero et al. 2002; Wilson et al. 2003; Krips 
et al. 2008). Krips et al. (2008) suggest that the gas 
densities in AGN host galaxies (< 10^ ® cm“®) are lower 
than starburst host galaxies ( 10 ® ° - 10 ® ® cm“®), and a 
common interpretation relies on a clumpy ISM near star¬ 
forming regions (which reduces the filling factor) and a 
continuous ISM near the AGN. Since our current obser¬ 
vations (^ 200 pc resolution) cover a significantly large 
area and the beam filling factor may be small (< 1 ’a at 
EO, El, and E2 are < 0.03, 0.04 - 0.06 and 0.01 - 0.04, 
respectively), higher resolution observations (< 0'.'5) are 
required to confirm this scenario. 

In addition, our modeling shows higher [HGNJ/pGO+j 
near the eastern nucleus (> 5) than that in the surround¬ 
ing clumps (< 4) and the overlap region (1-9). The el¬ 
evated [HCN]/[HCO'®] is explained by two mechanisms 
(Krips et al. 2008). One is far-UV radiation from OB 
stars in young starbursts (Sternberg & Dalgarno 1995), 
and the other is strong X-ray radiation from an AGN 
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Fig. 11.— (a) The Kennicutt-Scimidt law of VV 114 overlaid on other galaxies. Filled squares show regions of VV 114, while crosses 
show galaxies in the sample of Kennicutt (1998). The dashed line and the dotted line indicate the “starburst” sequence and “normal disk” 
sequence, respectively (Daddi et al. 2010). (b) Distribution of the box-averaged gas depletion time (rgas = Dh 2 /^SFr)- The gas depletion 
time in color scale ranges from 0-1 Gyr. The black crosses show the positions of the nuclei defined by the peak positions of the Ks-band 
observation (Tateuchi et al. 2012). The open square is a region which we cannot estimate the gas depletion time because of non-detection 
of the Paa emission (Tateuchi et al. 2012). (c) The Keniccutt-Schmidt law with the ^3_2/i-o- The ratio in color scale ranges from 0 to 
0.8. (d) The Keniccutt-Schmidt law with the -Ri 2 /i 3 - The ratio in color scale ranges from 0 to 40. 


(Maloney et al. 1996). Because of different penetrat¬ 
ing lengths between far-UV and X-ray emission, pho¬ 
ton dominated regions (PDRs) are created at the sur¬ 
face of gas clouds and X-rays penetrate deeply into the 
circumnuclear disk (CND), forming large X-ray domi¬ 
nated regions. As a consequence of this volume versus 
surface effect, the X-ray radiation from an AGN may 
produce higher HCN abundances than the UV radia¬ 
tion of starburst activities (Krips et al. 2008). To some 
degree, ionization effects from cosmic rays (Wild et al. 
1992) such as supernovae or strong shocks are suspected 
to significantly increase the HCO’*' abundance while po¬ 
tentially decreasing the HCN abundance, thus yielding 
lower i?HCN/HCO+ in evolved starbursts than in AGNs. 
The high [HCN]/[HCO“''] near the eastern nucleus and 
low [HCNJ/pCO+J and strong/extended 8 GHz contin¬ 
uum detection at the surrounding clumps (Condon et al. 
1991) are all consistent with a presence of an AGN in 
the eastern nucleus, surrounded by star-forming dense 
clumps. 


6 . 2 . Spatially Resolved Kennicutt-Sehmidt Law 


Observational studies of galaxies at global scales have 
shown that the surface density of SFR and that of cold 
gas traced in CO (1-0) obey a power law relation (KS 
law; Schmidt 1959; Kennicutt 1998). ULIRGs are sys- 
temically shifted from the normal galaxy population in 
the Esfr - Sh 2 phase (Komugi et al. 2005; Daddi et al. 
2010; Genzel et al. 2010; Leroy et al. 2013). It is sug¬ 
gested that systems lower in IR luminosity (e.g., LIRGs) 
occupy the region between the “starburst” sequence and 
the “normal disk” sequence in the KS law. Galaxies 
in the “starburst” sequence have shorter gas depletion 
time (xgas = Eh 2 /^SFr ^ 0.1 Gyr) relative to galaxies 
in the “normal disk” sequence (rgas ~ 1 Gyr; Daddi et al. 
2010; Bournaud et al. 2011). The spatially resolved sur¬ 
face densities of the SFR and the molecular gas mass of 
VV 114 are shown in Table 10 and Figure 11. The star¬ 
forming regions of VV 114 fill the gap between the “nor¬ 
mal disk” and “starburst” sequences (Figure 11a). We 
also show the spatial distribution of Tgas in Figure 11b. 
The data points close to the “starburst” sequence are lo¬ 
cated along the eastern nucleus (<0.2 ( 7 ^) Gyr) and 
the overlap region (= 0.2 - 0.4 (q^) Gyr), while those 
near the “normal disk” sequence are located in the NW 
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and SE arms (> 0.8 Gyr). The spatial distribution 

of SsFR and Eh 2 are consistent with the distributions of 
previous optical, UV, and X-ray studies (Alonso-Herrero 
et al. 2002; Le Floc’h et al. 2002; Grimes et al. 2006). 
Regions with higher EgpR and Eh 2 clearly show higher 
i? 3 _ 2 /i-o and i?i 2 /i 3 (Figures 11c and lid). 

In summary, transition from the “normal disk” to 
“starburst” sequence may occur when the molecular 
clouds become excited and dense at the nuclei and the 
overlap region. Moreover, gas clouds with high i? 3 _ 2 /i_o 
have high Esfr “ Sh 2 , and this is consistent with past 
studies which suggest that the i? 3 _ 2 /i_o correlates with 
the local Ha flux (Minamidani et al. 2008; Fuji! et al. 
2014). The R 12/13 also shows a similar trend, and this 
is also consistent with the past studies (> 20 in central 
kpc regions of U/LIRGs, 10 - 15 in normal starburst 
galaxies, and ^ 5 in Galactic GMCs; Aalto et al. 1997): 
The reason for the elevated i?i 2 /i 3 in starburst regions 
of VV 114 will be discussed in detail in $6.3. 

6.3. CO Isotope Ratio Enhancement in the Molecular 
“Filament” 

We suggest from our RADEX modelings that the east¬ 
ern nucleus and the overlap region have extremely high 
[CO]/[^^GO] (> 200), which is at least two times higher 
than the Galactic value (~ 70; Wilson & Rood 1994). 
The Paa peaks roughly coincide with the regions where 
high [CO]/[^^GO] are expected, suggesting that the in¬ 
creased [CO]/[^^CO] is related to the star formation ac¬ 
tivity. Similarly high values are seen in the overlap region 
of NGC4038/9 (Wilson et al. 2003) and the Taffy (Zhu 
et al. 2007). Zhu et al. (2007) suggested that the ex¬ 
treme [CO]/[^^CO] value in the bridge is explained by 
three scenarios, 1 ) selective isotope photodissociation in 
the diffuse clouds and shocked region, 2) CO enrichment 
around starburst activities, and/or 3) the destruction 
and recombination of molecules after shock. We briefly 
explain each scenario below, but our current data is in¬ 
sufficient for us to identify the exact cause of the high 
[CO]/pCO] in VV 114. 

The first possibility of [CO]/[^^CO] enhancement is 
the deficiency in ^^CO. Sheffer et al. (1992) suggest 
that selective isotope photodissociation can reduce the 
abundance in diffuse clouds, because CO is self- 
shielded to a greater extent. Thus, the ISM surrounding 
young starbursts and/or shocked regions show elevated 
[CO]/[^^CO] (Zhu et al. 2007). The ISM in the nuclei 
and the overlap region of VV 114 show extremely high 
[CO]/[^^CO], presumably due to intense starburst activ¬ 
ities and/or large-scale shocks. 

The second possibility is that massive stars end their 
life as supernovae and expel a large amount of in 
the interstellar medium. While the elemental abun¬ 
dances (e.g. C and S) are not directly related to the 
molecular abundances (e.g., CS; Casoli et al. 1992), once 
the synthesized elements are dispersed in the interstellar 
medium, molecules (e.g., CO, CS, and CN; Henkel et al. 
2014) form as soon as the temperature and density con¬ 
ditions are favorable. This occurs with a timescale of a 
few 10® yr (Langer & Graedel 1989). 

For the overlap region, the destruction and recombi¬ 
nation of molecules after shocks (see §6.5) are possible 
mechanisms to enhance the [CO]/[®®CO] (the third pos- 



Fig. 12.— Distribution of the box-averaged gas-to-dust ratio 
map. The black crosses show the positions of the nuclei defined 
by the peak positions of the Ks-band observation (Tateuchi et al. 
2012). The open squares are regions which we cannot estimate the 
gas-to-dust ratio because of non-detection of the 340 GHz contin¬ 
uum emission. 

sibility). The recombination timescale of H 2 and CO 
molecules after shock destruction are shorter than that 
of ^®CO, since ionized photons from shocked regions lead 
to selective isotope photodissociation (Zhu et al. 2003). 
Shielded regions from the radiation held are needed to 
form rare ^®CO (Abundant CO can form self-shielded 
regions). Moreover, the rare isotope molecules generally 
need a longer time to form, because collisions between 
molecules and dust grains are less frequent (Zhu et al. 
2007). 

6.4. Gas-to-Dust Ratio, M^i^/M^ust 

The gas-to-dust ratio, Mn^/M^ust, provides an im¬ 
portant measure of the relative abundance between gas 
and metallicity. The average M-a^/Mdust over the en¬ 
tire galaxy is often derived in single-dish work, and typ¬ 
ical M}d^/Mdust is 200 - 300 for local U/LIRGs (Contini 
& Contini 2003; Yao et al. 2003; Seaquist et al. 2004), 
and 15 - 231 in high-z sources (Solomon & Vanden Bout 
2005). Wilson et al. (2008) found MHa/Afdust = 357 ± 
95 from a sample of 13 U/LIRGs, including VV 114, ob¬ 
served at kpc resolution. 

We use the gas and dust masses derived in §5.4 to inves¬ 
tigate the distribution of Mnj/A/dust (Figure 12). The 

smallest value of (128 ± 16) (y^) (;^^) occurs in the 

eastern nucleus, which is similar to the Galactic value 
(100; Hildebrand 1983), while higher values of (371 ± 

118) (tf) and (339 ± 60) (ff) (^) occur 

in the western nucleus and the overlap region, respec¬ 
tively. The clear differences between the two nuclei may 
suggest a local gradient in the metallicity. For the over¬ 
lap region, cold dust associated with diffuse gas clouds 
cannot avoid the collision. This tends to increase the 
Mna/Vldust, because shocks destruct dust particles pref¬ 
erentially (Zhu et al. 2007). On the other hand, the low 
Mhj /Vfdust in the eastern nucleus may be due to intense 
starbursts producing dust-rich environments. 

6.5. Fractional Abundances of CS, CH 3 OH, and CN 

Table 15 shows the properties of the detected molecu¬ 
lar lines which are not used in the RADEX calculations. 
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Fig. 13.— Logarithmic fractional abundances relative to H 2 
([X]/[H 2 ]) of selected extragalactic sources compared to those of 
specific regions of VV 114, as presented in Table 15. The red, 
green, and blue columns show R18, R21a, and R39a, respectively. 
The pink, light blue, and yellow columns show NGC 253, M 82, 
and IC 342, respectively. Arrows represent upper limits. 

Either the dense gas component of VV 114 has extreme 
variations in excitation among the molecular clumps in 
the filament (see § 6 . 1 ), or there is widespread chemical 
differentiation across the filament. The fractional abun¬ 
dances [Vx/Vhj] of the different astrochemical species 
provide evidence of varying chemical influences due to 
star formation, physical conditions, and dynamics across 
the galaxy disks. We use the H 2 column densities, de¬ 
rived from the RADEX calculations, which are 10^°®, 
10^^’^, and 10^^’^ cm“^ at R18 (AGN), R21a (starburst), 
and R39a (overlap region), respectively. Column densi¬ 
ties of each molecules are determined by equation (7) 
assuming an optically thin emission under LTE. The Tex 
values determined from equation (3) are 38.7 ± 1.9 K, 
58.8 ± 2.9 K, and 52.6 ± 2.6 K at R18, R21a, and 
R39a, respectively. The derived [Vx/Vhj] are listed in 
Table 15. 

In Figure 13, we show the fractional abundances for 
CS, CH 3 OH, and CN in VV 114, and the same ratios for 
a sample of nearly galaxies, NGC 253, M82, and IC 342, 
taken from line surveys available in the literature (Henkel 
et al. 1988; Mauersberger et al. 1989; Huettemeister et al. 
1997; Martin et al. 2006). M82 has a relatively old star- 
burst at its core, with an average stellar population age 
of ~ 10 Myr (Konstantopoulos et al. 2009). This creates 
strong UV fields, therefore the PDR dominates its chem¬ 
istry (Aladro et al. 2011). Figure 13 shows that R39a 
has higher CH 3 OH abundance than M82, and small CS 
and CN abundances. A pure PDR similar to M82 may 
explain the molecular abundances we observe in R18. 

The molecular abundances for the overlap region and 
NGC 253 share similar characteristics. NGC 253 is 
thought to be in an early stage of starburst evolution, and 
has young stellar populations in its nucleus (~ 6 Myr; 
Fernandez-Ontiveros et al. 2009). The chemistry in the 
nucleus of NGC 253 is dominated by large-scale shocks 
(Aladro et al. 2011), and we suggest that the overlap re¬ 
gion of VV 114 is also dominated by shocks. The low 
i?HCN/HCO+ 9-1 the overlap region are further evidences 
for a shock dominated region (Krips et al. 2008). 

6 . 6 . Merger-driven Tidal Dwarf Galaxy Formation 


Tidal dwarf galaxies (TDGs) are gas-rich irregular 
galaxies made out of stellar and gaseous material pulled 
out by tidal forces from the disks of the colliding parent 
galaxies into the intergalactic medium. They are found 
at the ends of long tails and host active star-forming re¬ 
gions (Braine et al. 2000). Hibbard et al. (2001) and Gao 
et al. (2001) found the HI gas mass of 4.1 x 10® Mq and 
the molecular gas mass of 4 x 10® Mq at the edge of the 
southern tail of NGC 4038/9. 

We found an elevated i? 3 - 2 /i-o (0.36 ± 0.01), SFR 
(0.10 ± 0.05 Mq yr-i), and Mh^ (- 3.8 x 10^ {^) 
Mq) at the edge of the southern tidal arm (R38). The 
derived SFR and Mh 2 of R38 are comparable to those of 
TDG candidates in other galaxies (Braine et al. 2001). 
The gas depletion time of (0.40 ± 0.22) Gyr is 

shorter than the rest of the gas in the tidal arm (> 
0.5 (q^) Gyr). According to the RADEX modeling, 
while the ranges of Tkin and uhs are not confined well, 
the best fitting values (35 K, 10^ ® cm“®) are slightly 
higher than those in the middle of the tidal arm, R36a 
and R37a (25 ~ 30 K, 10^ ° - 10^'^ cm“®). We suggest 
that R38 is a forming tidal dwarf galaxy at the edge of 
the tidal arm of VV 114. Future high sensitivity optical 
and high resolution HI observations will allow us to con¬ 
strain the star formation and the atomic gas properties 
of R38. 

7. CONCLUSION 

We investigate the physical conditions of the molec¬ 
ular gas in the mid-stage merger VV 114. We present 
high-resolution observations of molecular gas and dust 
continuum emission in this galaxy using ALMA band 3 
and band 7. This study includes the first detection of ex- 
tranuclear CH 3 OH (2-1) emission in interacting galaxies. 
The results can be summarized as follows: 

1. We find that the CO (1-0) and CO (3-2) lines show 
significantly extended structures (i.e., the north¬ 
ern and southern tidal arms), the central filament 
across the galaxy disks, and double-peaks in the 
overlap region, while the ^®CO ( 1 - 0 ) line is only 
detected at the central filament. The filament 
is also identified by the strong CN (I 3/2 - O 1 / 2 ), 
HCO+ (4-3), 110 GHz, and 340 GHz continuum 
emission. 

2. Higher i? 3 _ 2 /i_o (0.5 — 0.8) and R 12/13 (20 — 50) 
are detected at the central filament. These higher 
ratios indicate that the central filament has highly 
excited (but not thermalized) molecular ISM, and 
the eastern nucleus is nearly thermalized when it 
is observed with a 800 pc beam. 

3. The unresolved eastern nucleus has the highest 
i?HCN/HCO+ (1-34 ± 0.09), while the dense gas 
clumps near the eastern nucleus have significantly 
lower values (~ 0.5). The broad HCN (4-3) and 
HCO+ (4-3) (^ 290 km s“^) emission lines seen 
in the unresolved eastern nucleus suggests an ob¬ 
scured AGN (see also paper I). 

4. Radiative transfer analysis of the CO (1-0), CO (3- 
2 ), and ^®CO ( 1 - 0 ) emission enables us to map 
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physical parameters of the “diffuse” gas of an in¬ 
teracting LIRG with 800 pc scale for the first 
time. The analysis suggests that “diffuse” gas 
clouds in the filament have warmer/denser condi¬ 
tions than those in the galaxy disks. This is con¬ 
sistent with predictions from merger simulations. 
Our analysis also suggest that the [CO]/[^^CO] is 
enhanced in the central filament. The extremely 
high [CO]/[^^CO] values are more important than 
the moderately optically thick ( 1 - 0 ) emission 

to explain the high R 12/13 in VV 114. 

5. Radiative transfer analysis of the HCN (4-3), 

HCO+ (4-3), and (3-2) allow us to com¬ 

pare the dense gas clouds around AGN, starburst 
activities, and the overlap region. These results 
show that dense gas clouds around AGN have tt-Hj 
= 10'5 ° - lO'^-^ cm-3 and Tkin > 100 K with 
[HCN]/[HCO+] > 5, while gas clumps around star- 
burst activities show riHa = 10 ^ ® - 10 ^ ® cm“^ and 
Tkin = 40 -100 K with [HCN]/[HCO+] < 4. In ad¬ 
dition, the analysis shows that the overlap region 
has nH 2 = 10® ° - 10® ° cm'^ and = 5 - 90 K 
with [HGN]/[HGO+] = 1-9. 

6 . The spatially resolved Kennicutt-Schmidt law in 
VV 114 clearly connects the “starburst” sequence 
with the “normal disk” sequence. Most of the data 
points near the “starburst” sequence are found in 
the nuclei and the overlap region, whereas the data 
points near the “normal disk” sequence are found 
in the tidal arms. We also find the i? 3 _ 2 /i-o and 
i?i 2 /i 3 are well correlated with the Ssfr- 

7. The A/H,/Mdu.t of (128 ± 16) (ff) (^) in 

the eastern nucleus of VV 114 is comparable to 
the Galactic value, but it is a factor of two 
higher than that in the overlap region of (339 ± 

60) (q^) Since the 340 GHz emission is 

spatially correlated with dense gas tracers, the cold 
dust in VV 114 appears to be closely related to the 
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TABLE 1 

Log of alma Observations 


UT date 

( 1 ) 

Spectral windows 

LSB USB 

[GHz] [GHz] 

(2) (3) 

Configuration 

.Nant Array Auaseline 
[m] 

(4) (5) ( 6 ) 

Tsys 

[K] 

(7) 

MRS Amplitude caibrator 

[arcsec.] 

( 8 ) (9) 

robs 

[min.] 

( 10 ) 

2011 Nov 6 

101.5, 103.5 

114.0, 115.1 

16 

CMP 

18 - 196 

65 - 89 

18 

Uranus 

41 

2012 May 4 

101.5, 103.5 

114.0, 115.1 

15 

EXT 

39 - 402 

48 - 62 

8 

Neptune 

40 

2012 Mar 27 

97.5, 99.5 

110.2, 111.5 

17 

EXT 

18 - 401 

54 - 73 

19 

Neptune 

22 

2012 Jul 2 

97.5, 99.5 

110.2, 111.5 

20 

EXT 

16 - 402 

71 - 117 

21 

Neptune 

39 

2012 Nov 5 

331.1, 333.0 

343.5, 345.3 

14 

CMP 

12 - 135 

125 - 172 

9 

Uranus 

66 

2012 Nov 5 

331.1, 333.0 

343.5, 345.3 

14 

CMP 

12 - 135 

108 - 155 

9 

Uranus 

67 

2012 Nov 5 

331.1, 333.0 

343.5, 345.3 

14 

CMP 

12 - 135 

124 - 175 

9 

Callisto 

67 

2012 Jun 1 

342.0, 344.0 

354.5, 356.0 

18 

EXT 

15 - 402 

150 - 213 

7 

Uranus 

78 

2012 Jun 2 

342.0, 344.0 

354.5, 356.0 

20 

EXT 

15 - 402 

108 - 160 

7 

Uranus 

80 

2012 Jun 3 

342.0, 344.0 

354.5, 356.0 

20 

EXT 

15 - 402 

103 - 130 

7 

Uranus 

45 

Note. — Column 2 and 3: 

Central frequencies of the spectral windows (spw). All spw have the frequency coverage of 1.875 GHz. 


Column 4: Number of available antennas. Column 5: ALMA antenna configuration. CMP is the compact configuration and EXT is 
the extended configuration. Column 6: Range of projected length of baselines for VV 114. Column 7: DSB system temperature toward 
VV 114. Column 8: Maximum recoverable scale (MRS) of the configuration. This is defined by ~ 0.6 A/(minimum Lbaseline)- Column 
9: Observed calibrators for amplitude correction. Column 10: Total integration time on the galaxy. 


TABLE 2 

ALMA Observational Properties 


Emission 

( 1 ) 

Band 

( 2 ) 

t^rest 

[GHz] 

( 3 ) 

Beam size 
[arcsecond] 

( 4 ) 

P.A. 

[deg] 

( 5 ) 

AV 

[km s“l] 
( 6 ) 

Noise rms 
[mjy bearn”^] 

( 7 ) 

[mK] 

( 8 ) 

CH 3 OH {2k-lk) 

3 

96.74 

2.03 

X 

1.34 

85.7 

30 

1.0 

46 

CS (2-1) 

3 

97.98 

2.01 

X 

1.37 

83.6 

30 

0.9 

40 

i^CO ( 1 - 0 ) 

3 

110.20 

1.77 

X 

1.20 

85.8 

30 

1.0 

46 

CN (I1/2-O1/2) 

3 

113.14 

1.97 

X 

1.27 

-85.8 

30 

1.0 

37 

CN (I3/2-O1/2) 

3 

113.49 

1.98 

X 

1.29 

-84.7 

30 

1.1 

39 

CO (1-0) 

3 

115.27 

1.97 

X 

1.35 

82.3 

10 

2.3 

76 

CS (7-6) 

7 

342.88 

0.47 

X 

0.39 

54.2 

30 

0.7 

38 

CO (3-2) 

7 

345.80 

1.64 

X 

1.17 

112.6 

10 

2.1 

11 

HCN (4-3) 

7 

354.51 

0.46 

X 

0.38 

51.5 

30 

0.8 

42 

HCO+ (4-3) 

7 

356.73 

0.45 

X 

0.37 

53.4 

30 

0.9 

50 

Continuum 

3 

110 

1.89 

X 

1.28 

81.8 


0.05 

2.1 

Continuum 

7 

340 

1.33 

X 

1.12 

119.6 


0.11 

0.8 

Continuum 

7 

340 

0.45 

X 

0.38 

56.2 


0.07 

4.3 


Note. — Column 1: Identified emission. Column 2: Band which includes the molecular line and 
continuum emission. Column 3: Rest frequency of the line or mean frequency of the continuum. 
Column 4: Major and minor axes (FWHM) of the synthesized beam. Column 5: Position angle of 
the synthesized beam. Column 6: Velocity resolution of our binning images. Column 7 and 8: Noise 
rms intensity in the data which have velocity resolutions shown in Column 6. The noise in Column 
8 is in Rayleigh-Jeans brightness temperature. 
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TABLE 3 

Radio Continuum Flux of VV 114 


ID 

(1) 


00 


Siio 

[mjy] 

(3) 


5340 

[mJy] 

(4) 

R7 

0.62 

± 

0.08 


< 

0.18 

0.58 

± 

0.14 

R8 

0.70 

± 

0.08 


< 

0.18 


< 

0.42 

R9 

0.82 

± 

0.08 

0.21 

± 

0.06 


< 

0.42 

RIO 

0.81 

± 

0.08 

0.20 

± 

0.06 


< 

0.42 

Rll 

0.92 

± 

0.08 


< 

0.18 


< 

0.42 

R17 

0.42 

± 

0.08 


< 

0.18 


< 

0.42 

R18t 

3.37 

± 

0.08 

1.48 

± 

0.06 

5.17 

± 

0.14 

R19 

2.40 

± 

0.08 

1.03 

± 

0.06 

3.03 

± 

0.14 

R20 

1.14 

± 

0.08 

0.54 

± 

0.06 

2.41 

± 

0.14 

R2ltt 

5.00 

± 

0.08 

1.86 

± 

0.06 

8.20 

± 

0.14 

R22 

1.85 

± 

0.08 

0.81 

± 

0.06 

4.37 

± 

0.14 

R23 

1.34 

± 

0.08 

0.50 

± 

0.06 

3.55 

± 

0.14 

R24 

0.67 

± 

0.08 


< 

0.18 

1.19 

± 

0.14 

R25 

1.07 

± 

0.08 

0.39 

± 

0.06 

3.01 

± 

0.14 

R26 

0.81 

± 

0.08 

0.33 

± 

0.06 

1.73 

± 

0.14 

R27 

1.29 

± 

0.08 

0.41 

± 

0.06 

1.40 

± 

0.14 

R28 

0.81 

± 

0.08 

0.21 

± 

0.06 

0.76 

± 

0.14 

R29 

0.35 

± 

0.08 


< 

0.18 

0.71 

± 

0.14 

R30 


< 

0.24 


< 

0.18 

0.70 

± 

0.14 

R39ttt 

1.18 

± 

0.08 

0.21 

± 

0.06 

3.35 

± 

0.14 


Note. —Column 2: 8.44 GHz continuum flux (Con¬ 
don et al. 1991). Column 3: 110 GHz continuum 

flux obtained by ALMA/band 3. Column 4: 340 GHz 
continuum flux obtained by ALMA/band 7.; We only 
show the statistical error in this table. The system¬ 
atic error of absolute flux calibration is estimated to 
be ~ 5% in band 3 and ~ 10% in band 7. frepresents 
boxes contained the obscured AGN deflned by paper 1. 
tfrepi’esents boxes contained the nuclear starbursts de¬ 
fined by paper 1. ftfrepresents boxes at the overlap 
region. 
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TABLE 4 

CO (1-0), (1-0), AND (3-2) Integrated Intensities 


ID 

12 C 0 ( 1 - 0 ) 

1=00 (3-2) 

13 CO (1-0) 

^3 

-2 

/l-O 

^12/13 


[Jy km 

s-l] 

[Jy km 

s-l] 

[Jy km 

s-l] 







( 1 ) 



( 2 ) 



(3) 



(4) 



(5) 


( 6 ) 

R1 

2.12 

± 

0.06 

5.92 

± 

0.08 




0.31 

± 

0.01 




R2 

1.41 

± 

0.05 

7.47 

± 

0.06 

0.29 

± 

0.02 

0.59 

± 

0.02 

4 

± 

1 

R3 

3.62 

± 

0.09 

5.45 

± 

0.09 




0.17 

± 

0.01 




R4 

3.54 

± 

0.08 

8.69 

± 

0.10 

0.25 

± 

0.04 

0.27 

± 

0.01 

13 

± 

2 

R5 

3.46 

± 

0.08 

9.16 

± 

0.09 

0.40 

± 

0.03 

0.29 

± 

0.01 

7 

± 

1 

R 6 

3.66 

± 

0.10 

10.71 

± 

0.11 




0.33 

± 

0.01 




R7 

5.26 

± 

0.10 

19.64 

± 

0.12 

0.40 

± 

0.03 

0.41 

± 

0.01 

12 

± 

1 

R 8 

3.71 

± 

0.09 

13.44 

± 

0.12 

0.38 

± 

0.03 

0.40 

± 

0.01 

8 

± 

1 

R9 

9.59 

± 

0.11 

33.88 

± 

0.12 

0.43 

± 

0.04 

0.39 

± 

0.01 

20 

± 

2 

RIO 

12.69 

± 

0.11 

47.74 

± 

0.12 

0.37 

± 

0.05 

0.42 

± 

0.01 

32 

± 

4 

Rll 

16.26 

± 

0.13 

58.93 

± 

0.14 

0.49 

± 

0.04 

0.40 

± 

0.01 

30 

± 

3 

R12 

4.61 

± 

0.10 

13.53 

± 

0.10 




0.33 

± 

0.01 




R13 

5.22 

± 

0.10 

16.02 

± 

0.11 

0.38 

± 

0.05 

0.34 

± 

0.01 

12 

± 

2 

R14 

4.70 

± 

0.09 

12.22 

± 

0.09 

0.35 

± 

0.03 

0.29 

± 

0.01 

12 

± 

1 

R15 

5.03 

± 

0.09 

10.62 

± 

0.09 

0.35 

± 

0.03 

0.23 

± 

0.01 

13 

± 

1 

R16 

5.59 

± 

0.09 

11.23 

± 

0.08 

0.38 

± 

0.02 

0.22 

± 

0.01 

14 

± 

1 

R17 

6.44 

± 

0.10 

14.33 

± 

0.11 

0.38 

± 

0.05 

0.25 

± 

0.01 

15 

± 

2 

R18t 

16.54 

± 

0.13 

106.39 

± 

0.15 

0.60 

± 

0.05 

0.71 

± 

0.01 

25 

± 

2 

R19 

11.09 

± 

0.12 

54.98 

± 

0.14 

0.60 

± 

0.05 

0.55 

± 

0.01 

17 

± 

2 

R20 

18.33 

± 

0.13 

93.14 

± 

0.14 

0.64 

± 

0.06 

0.56 

± 

0.01 

26 

± 

2 

R2ltt 

33.04 

± 

0.15 

225.47 

± 

0.16 

1.09 

± 

0.07 

0.76 

± 

0.01 

28 

± 

2 

R22 

28.67 

± 

0.15 

162.41 

± 

0.16 

0.90 

± 

0.07 

0.63 

± 

0.01 

29 

± 

2 

R23 

33.65 

± 

0.16 

177.02 

± 

0.17 

0.95 

± 

0.07 

0.58 

± 

0.01 

32 

± 

2 

R24 

16.12 

± 

0.14 

68.28 

± 

0.15 

0.48 

± 

0.07 

0.47 

± 

0.01 

31 

± 

4 

R25 

25.23 

± 

0.16 

138.20 

± 

0.16 

0.86 

± 

0.06 

0.61 

± 

0.01 

27 

± 

2 

R26 

20.44 

± 

0.14 

97.56 

± 

0.15 

0.69 

± 

0.06 

0.53 

± 

0.01 

27 

± 

3 

R27 

12.45 

± 

0.12 

67.10 

± 

0.13 

0.54 

± 

0.06 

0.60 

± 

0.01 

21 

± 

2 

R28 

5.78 

± 

0.10 

36.14 

± 

0.14 

0.35 

± 

0.03 

0.69 

± 

0.01 

15 

± 

1 

R29 
















R30 

6.39 

± 

0.12 

30.50 

± 

0.12 

0.28 

± 

0.03 

0.53 

± 

0.01 

21 

± 

3 

R31 

4.93 

± 

0.11 

19.03 

± 

0.13 

0.30 

± 

0.04 

0.43 

± 

0.01 

15 

± 

2 

R32 

2.34 

± 

0.05 

7.34 

± 

0.08 

0.21 

± 

0.02 

0.35 

± 

0.01 

10 

± 

1 

R33 

1.44 

± 

0.05 

2.03 

± 

0.05 




0.16 

± 

0.01 




R34 

2.81 

± 

0.08 

11.65 

± 

0.10 

0.32 

± 

0.03 

0.46 

± 

0.01 

8 

± 

1 

R35 

5.36 

± 

0.09 

9.40 

± 

0.09 

0.30 

± 

0.03 

0.19 

± 

0.01 

16 

± 

2 

R36 

5.55 

± 

0.12 

17.95 

± 

0.14 

0.31 

± 

0.02 

0.36 

± 

0.01 

17 

± 

1 

R37 

4.84 

± 

0.11 

8.65 

± 

0.11 

0.31 

± 

0.03 

0.20 

± 

0.01 

14 

± 

2 

R38 

2.66 

± 

0.07 

8.51 

± 

0.09 

0.33 

± 

0.02 

0.36 

± 

0.01 

7 

± 

1 

R39ttt 

27.54 

± 

0.15 

155.31 

± 

0.15 

0.97 

± 

0.06 

0.63 

± 

0.01 

26 

± 

2 


Note. — Column 1: These numbers are labeled at the ratio map of Figure 7. 
Column 2: Integrated ^^CO (1—0) intensity at an emission region. Column 3: Inte¬ 
grated ^^CO (3—2) intensity at an emission region. Column 4: Integrated ^^CO (3—2) 
intensity at an emission region. Column 5: The ^^CO (3—2)/CO (1—0) integrated 
intensity ratio. Column 6: The ^^CO (1—0)/^^CO (1—0) integrated intensity ratio.; 
We only show the statistical error in this table. The systematic error of absolute 
flux calibration is estimated to be ~ 5% in band 3 and ~ 10% in band 7. frepresents 
boxes contained the obscured AGN defined by paper I. ffrepresents boxes contained 
the nuclear starbursts defined by paper I. ftfrepresents boxes at the overlap region. 
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TABLE 5 

HCN (4-3) AND HCO+ (4-3) Integrated 
Intensities 


ID HCN (4-3) HCO+ (4-3) Rhcn/hco+ 

[Jy km s“^] [Jy km s“^] 

(1) (2) (3) (4) 


sot 

1.34 

± 

0.09 

0 .88± 

0.10 

1.52 

± 

0.20 

SI 

0.11 

± 

0.03 

0.34± 

0.06 

0.34 

± 

0.10 

S2 


< 

0.05 

0 .12± 

0.06 


< 

0.43 

S3tt 

0.49 

± 

0.07 

1.05± 

0.08 

0.46 

± 

0.07 

S4 

0.13 

± 

0.06 

0.63± 

0.08 

0.20 

± 

0.09 

S5tt 

0.81 

± 

0.07 

2.36± 

0.09 

0.34 

± 

0.03 

S6 

0.40 

± 

0.07 

0.92± 

0.08 

0.43 

± 

0.09 

S7 

0.19 

± 

0.07 

0.42± 

0.08 

0.45 

± 

0.18 

S8 

0.06 

± 

0.04 

0.14± 

0.06 

0.44 

± 

0.31 

S9 

0.07 

± 

0.05 

0.13± 

0.06 

0.54 

± 

0.46 

SIO 


< 

0.04 

0.28± 

0.07 


< 

0.14 

Sllttt 


< 

0.05 

0.15± 

0.06 


< 

0.33 

S12ttt 


< 

0.04 

0.32± 

0.07 


< 

0.12 

S13 


< 

0.04 

0.17± 

0.06 


< 

0.24 

S14 


< 

0.05 

0.19± 

0.07 


< 

0.26 


Note. — Column 1: These numbers are labeled at the 
ratio map of Figure 7. Column 2: Integrated HCN (4— 
3) intensity at an emission region. Column 3: Integrated 
HCO'*’ (4—3) intensity at an emission region. Column 4: 
The HCN (4—3)/HCO"^ (4—3) integrated intensity ratio.; 
We only show the statistical error in this table. The sys¬ 
tematic error of absolute flux calibration is estimated to 
be ~ 5% in band 3 and ~ 10% in band 7. frepresents 
boxes contained the obscured AGN defined by paper I. 
tfrepresents boxes contained the nuclear starbursts de¬ 
fined by paper I. fffrepresents boxes at the overlap re¬ 
gion. 






22 


TOSHIKI SAITO et al. 


TABLE 6 

Peak Brightness Temperature oe (1-0), (1-0), and (3-2) Emission 


ID 

Peak Ti, co (i_o) 
[K] 

Peak Tb.CO ( 3 - 2 ) 
[K] 

Peak Tb 1300 ( 1 - 0 ) 
[K] 

Peak /? 3 _ 2 /i_o 

Peak -R 12/13 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


R 1 

0.40 

± 

0.07 

0.13 

± 

0.01 


< 

0.09 

0.32 

± 

0.07 

> 

4 

R2 

0.69 

± 

0.07 

0.25 

± 

0.01 


< 

0.09 

0.37 

± 

0.04 

> 

8 

R3 

0.61 

± 

0.07 

0.07 

± 

0.01 


< 

0.09 

0.12 

± 

0.02 

> 

7 

R4 

0.98 

± 

0.07 

0.30 

± 

0.01 


< 

0.09 

0.30 

± 

0.02 

> 11 

R5 

1.19 

± 

0.07 

0.42 

± 

0.01 

0.06 

± 

0.03 

0.35 

± 

0.02 

20 ± 

2 

R 6 

0.87 

± 

0.07 

0.21 

± 

0.01 


< 

0.09 

0.24 

± 

0.02 

> 10 

R7a 

1.34 

± 

0.07 

0.32 

± 

0.01 


< 

0.09 

0.24 

± 

0.01 

> 15 

R7b 

0.45 

± 

0.07 

0.36 

± 

0.01 

0.06 

± 

0.03 

0.79 

± 

0.12 

7 ± 

1 

R 8 a 

0.21 

± 

0.07 


< 

0.03 


< 

0.09 


> 

0.14 

> 

2 

R 8 b 

0.96 

± 

0.07 

0.24 

± 

0.01 

0.08 

± 

0.03 

0.25 

± 

0.02 

11 ± 

1 

R9a 

0.40 

± 

0.07 

0.12 

± 

0.01 


< 

0.09 

0.30 

± 

0.06 

> 

4 

R9b 

2.60 

± 

0.07 

0.99 

± 

0.01 

0.08 

± 

0.03 

0.38 

± 

0.01 

33 ± 

2 

RlOa 

0.80 

± 

0.07 

0.24 

± 

0.01 

0.05 

± 

0.03 

0.30 

± 

0.03 

17 ± 

2 

RlOb 

2.81 

± 

0.07 

1.27 

± 

0.01 

0.09 

± 

0.03 

0.45 

± 

0.01 

30 ± 

1 

Rlla 

2.11 

± 

0.07 

0.82 

± 

0.01 

0.12 

± 

0.03 

0.39 

± 

0.01 

18 ± 

1 

Rllb 

3.34 

± 

0.07 

1.36 

± 

0.01 

0.04 

± 

0.03 

0.41 

± 

0.01 

89 ± 

7 

R12a 

0.55 

± 

0.07 

0.23 

± 

0.01 


< 

0.09 

0.41 

± 

0.06 

> 

6 

R12b 

1.30 

± 

0.07 

0.41 

± 

0.01 


< 

0.09 

0.32 

± 

0.02 

> 14 

R13a 

1.12 

± 

0.07 

0.34 

± 

0.01 

0.06 

± 

0.03 

0.30 

± 

0.02 

18 ± 

1 

R13b 

0.67 

± 

0.07 

0.30 

± 

0.01 


< 

0.09 

0.45 

± 

0.05 

> 

7 

R14 

1.61 

± 

0.07 

0.39 

± 

0.01 

0.08 

± 

0.03 

0.24 

± 

0.01 

20 ± 

1 

R15 

1.77 

± 

0.07 

0.42 

± 

0.01 

0.08 

± 

0.03 

0.24 

± 

0.01 

22 ± 

1 

R16 

1.38 

± 

0.07 

0.33 

± 

0.01 

0.12 

± 

0.03 

0.24 

± 

0.01 

12 ± 

1 

R17 

1.64 

± 

0.07 

0.35 

± 

0.01 

0.05 

± 

0.03 

0.21 

± 

0.01 

36 ± 

3 

R18t 

3.52 

± 

0.07 

1.79 

± 

0.01 

0.14 

± 

0.03 

0.51 

± 

0.01 

25 ± 

1 

R19a 

1.01 

± 

0.07 

0.54 

± 

0.01 

0.09 

± 

0.03 

0.54 

± 

0.04 

11 ± 

1 

R19b 

1.99 

± 

0.07 

0.98 

± 

0.01 

0.09 

± 

0.03 

0.49 

± 

0.02 

23 ± 

1 

R20a 

3.33 

± 

0.07 

1.99 

± 

0.01 

0.11 

± 

0.03 

0.60 

± 

0.01 

30 ± 

1 

R20b 

1.13 

± 

0.07 

0.67 

± 

0.01 


< 

0.09 

0.59 

± 

0.04 

> 13 

R21att 

5.53 

± 

0.07 

4.49 

± 

0.01 

0.23 

± 

0.03 

0.81 

± 

0.01 

24 ± 

1 

R21btt 

1.22 

± 

0.07 

0.95 

± 

0.01 


< 

0.09 

0.78 

± 

0.05 

> 14 

R 21 ctt 

0.40 

± 

0.07 

0.72 

± 

0.01 


< 

0.09 

1.80 

± 

0.05 

> 

4 

R22a 

3.60 

± 

0.07 

1.93 

± 

0.01 

0.12 

± 

0.03 

0.54 

± 

0.01 

29 ± 

1 

R22b 

3.83 

± 

0.07 

2.61 

± 

0.01 

0.16 

± 

0.03 

0.68 

± 

0.01 

24 ± 

1 

R23a 

4.17 

± 

0.07 

2.56 

± 

0.01 

0.15 

± 

0.03 

0.61 

± 

0.01 

28 ± 

1 

R23b 

3.39 

± 

0.07 

1.87 

± 

0.01 

0.12 

± 

0.03 

0.55 

± 

0.01 

28 ± 

1 

R24a 

1.23 

± 

0.07 

0.46 

± 

0.01 

0.05 

± 

0.03 

0.37 

± 

0.02 

23 ± 

2 

R24b 

1.60 

± 

0.07 

0.96 

± 

0.01 

0.06 

± 

0.03 

0.60 

± 

0.03 

28 ± 

2 

R24c 

1.78 

± 

0.07 

0.58 

± 

0.01 

0.06 

± 

0.03 

0.32 

± 

0.01 

32 ± 

2 

R25a 

4.89 

± 

0.07 

3.18 

± 

0.01 

0.23 

± 

0.03 

0.65 

± 

0.01 

21 ± 

1 

R25b 

1.77 

± 

0.07 

0.59 

± 

0.01 

0.06 

± 

0.03 

0.33 

± 

0.01 

31 ± 

2 

R26a 

4.49 

± 

0.07 

2.57 

± 

0.01 

0.14 

± 

0.03 

0.57 

± 

0.01 

33 ± 

1 

R26b 

1.07 

± 

0.07 

0.31 

± 

0.01 

0.07 

± 

0.03 

0.29 

± 

0.02 

15 ± 

1 

R27a 

1.96 

± 

0.07 

1.31± 

0.01 

0.04 

± 

0.03 

0.67 

± 

0.02 

47 ± 

4 

R27b 

1.53 

± 

0.07 

0.75 

± 

0.01 

0.13 

± 

0.03 

0.49 

± 

0.02 

12 ± 

1 

R28a 

0.76 

± 

0.07 

0.37 

± 

0.01 

0.05 

± 

0.03 

0.48 

± 

0.05 

16 ± 

2 

R28b 

0.86 

± 

0.07 

0.54 

± 

0.01 


< 

0.09 

0.62 

± 

0.05 

> 10 

R29 


< 

0.21 


< 

0.03 


< 

0.09 






R30a 

0.77 

± 

0.07 

0.39 

± 

0.01 


< 

0.09 

0.51 

± 

0.05 

> 

9 

R30b 

1.39 

± 

0.07 

0.69 

± 

0.01 


< 

0.09 

0.50 

± 

0.03 

> 15 

R31a 

0.72 

± 

0.07 

0.33 

± 

0.01 

0.06 

± 

0.03 

0.46 

± 

0.03 

13 ± 

1 

R31b 

0.55 

± 

0.07 

0.24 

± 

0.01 


< 

0.09 

0.44 

± 

0.06 

> 

6 

R32 

0.55 

± 

0.07 

0.29 

± 

0.01 


< 

0.09 

0.52 

± 

0.07 

> 

6 

R33 

0.35 

± 

0.07 

0.14 

± 

0.01 


< 

0.09 

0.41 

± 

0.09 

> 

4 

R34 

0.83 

± 

0.07 

0.37 

± 

0.01 

0.05 

± 

0.03 

0.45 

± 

0.04 

18 ± 

2 

R35 

1.42 

± 

0.07 

0.30 

± 

0.01 


< 

0.09 

0.21 

± 

0.01 

> 16 

R36a 

0.65 

± 

0.07 

0.24 

± 

0.01 

0.07 

± 

0.03 

0.37 

± 

0.04 

9 ± 

1 

R36b 

0.47 

± 

0.07 

0.22 

± 

0.01 


< 

0.09 

0.47 

± 

0.07 

> 

5 

R37a 

0.79 

± 

0.07 

0.15 

± 

0.01 

0.06 

± 

0.03 

0.19 

± 

0.02 

14 ± 

1 

R37b 

0.71 

± 

0.07 

0.18 

± 

0.01 


< 

0.09 

0.26 

± 

0.03 

> 

8 

R38 

0.38 

± 

0.07 

0.16 

± 

0.01 

0.04 

± 

0.03 

0.42 

± 

0.08 

9 ± 

2 

R39attt 

4.91 

± 

0.07 

3.15 

± 

0.01 

0.24 

± 

0.03 

0.64 

± 

0.01 

21 ± 

1 

R39bttt 

1.96 

± 

0.07 

0.81 

± 

0.01 

0.06 

± 

0.03 

0.41 

± 

0.02 

31 ± 

2 


Note. — Column 1: These numbers are labeled at the ratio map of Figure 7. Column 2: Peak ^^CO (1—0) 
brightness temperature at an emission region. Column 3: Peak ^^CO (3—2) brightness temperature at an emission 
region. Column 4: Peak ^^CO (3—2) brightness temperature at an emission region. Column 5: The ^^CO (3— 
2)/CO (1—0) brightness temperature ratio. Column 6: The ^^CO (1—0)/^^CO (1—0) brightness temperature ratio.; 
We only show the statistical error in this table. The systematic error of absolute flux calibration is estimated to 
be ~ 5 % in band 3 and ~ 10% in band 7. trepresents boxes contained the obscured AGN deflned by paper 1. 
tfrepresents boxes contained the nuclear starbursts deflned by paper I. ftfrepresents boxes at the overlap region. 
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TABLE 7 


Peak Brightness Temperature of HCN (4-3) and 
HCO+ (4-3) 


ID Tb.HCN (4-3) 'Ib,HCO+(4-3) ^HCN/HCO + 

[K] [K] 

(1) (2) (3) (4) 


sot 

0.32 

± 

0.05 

0.27 

± 

0.05 

1.18 

± 

0.29 

SI 


< 

0.15 

0.20 

± 

0.05 


< 

0.75 

S2 


< 

0.15 

0.16 

± 

0.05 


< 

0.94 

S3tt 

0.20 

± 

0.05 

0.65 

± 

0.05 

0.31 

± 

0.08 

S4 

0.10 

± 

0.05 

0.51 

± 

0.05 

0.19 

± 

0.10 

S5tt 

0.32 

± 

0.05 

0.91 

± 

0.05 

0.35 

± 

0.06 

S6 

0.14 

± 

0.05 

0.55 

± 

0.05 

0.26 

± 

0.09 

S7 


< 

0.15 

0.24 

± 

0.05 


< 

0.63 

S8 


< 

0.15 

0.14 

± 

0.05 


< 

1.07 

S9 


< 

0.15 

0.17 

± 

0.05 


< 

0.88 

SIO 

0.10 

± 

0.05 

0.17 

± 

0.05 

0.60 

± 

0.34 

Sllttt 


< 

0.15 

0.18 

± 

0.05 


< 

0.83 

S12ttt 


< 

0.15 

0.24 

± 

0.05 


< 

0.63 

S13 


< 

0.15 

0.20 

± 

0.05 


< 

0.75 

S14 


< 

0.15 

0.23 

± 

0.05 


< 

0.65 


Note. — Column 1: These numbers are labeled at the 
ratio map of Figure 7. Column 2: Peak HCN (4—3) bright¬ 
ness temperature at an emission region. Column 3: Peak 
HCO'*" (4—3) brightness temperature at an emission region. 
Column 4: The HCN (4—3)/HCO'^ (4—3) brightness temper¬ 
ature ratio.; We only show the statistical error in this table. 
The systematic error of absolute flux calibration is estimated 
to be ~ 5 % in band 3 and ~ 10% in band 7. frepresents boxes 
contained the obscured AGN defined by paper I. ffrepresents 
boxes contained the nuclear starbursts defined by paper I. 
tffrepresents boxes at the overlap region. 


TABLE 8 

Peak Brightness Temperature of (3-2), HCN (4-3), and HCO+ (4-3) 


ID 

Peak Tt, CO (3-2) 
[K] 
(2) 

Peak Tij^hcn (4-3) 
[K] 
(3) 

Peak 7b jjco+(4-3) 
[K] 
(4) 

Peak 7 Jhco+/CO 

Peak 7 ?hcn/hco+ 

(1) 

(5) 

(6) 

EOt 

2.94 ± 0.04 

0.06 ± 0.01 

0.04 ± 0.01 

0.014 ± 0.004 

1.55 ± 0.49 

Eltt 

5.82 ± 0.04 

0.11 ± 0.01 

0.32 ± 0.01 

0.056 ± 0.002 

0.36 ± 0.03 

E2ttt 

4.71 ± 0.04 

0.04 ± 0.01 

0.08 ± 0.01 

0.017 ± 0.002 

0.52 ± 0.14 


Note. — Column 1: These numbers are regions where were convolved to the 1^2 X iVO resolution (P.A. = 119 deg.). 
Column 2: Peak ^^CO (3—2) brightness temperature at an emission region. Column 3: Peak HCN (4—3) brightness 
temperature at an emission region. Column 4: Peak HCO'*' (4—3) brightness temperature at an emission region. Column 
5: The HCO’*' (4—3)/CO (3—2) brightness temperature ratio. Column 6: The HCN (4—3)/HCO'*' (4—3) brightness 
temperature ratio.; We only show the statistical error in this table. The systematic error of absolute flux calibration is 
estimated to be ~ 5% in band 3 and ~ 10% in band 7. frepresents boxes contained the obscured AGN defined by paper I. 
ffrepresents boxes contained the nuclear starbursts defined by paper I. fffrepresents boxes at the overlap region. 


TABLE 9 

Parameters Used for RADEX Modelings 


Case 

Tkin 

log nn. 

log 7V(H2) 

abundance ratio 

box size 


[K] 

[cm-^] 

[cm-2l 


[pc] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

1 

5 - 300 (5) 

2 - 5 (0.1) 

18 - 22 (0.1) 

70 

800 

2 

5 - 400 (5) 

3 - 7 (0.1) 

21.2, 21.6, 21.5 

1 - 10 (1) 

320 


Note. — Column 2, 3, 4, 5: Fitting ranges (steps) of Tkini ^H2 > ^(^ 2 ): and 
abundance ratios. The abundance ratios mean the [CO]/[^^CO] and [HCNj/fHCO*^] 
in case 1 and 2, respectively. We fixed the [CO]/[ ^CO] in case 1 (Galactic value; 
Wilson &: Rood 1994) and the A^(H 2 ) (EO, El, and E2, respectively) in case 2. Column 
6: Box sizes.; We consider the statistical error and the systematic error in these 
calculations. The systematic error of absolute flux calibration is estimated to be ~ 
5% in band 3 and ~ 10% in band 7. 
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TABLE 10 

Gas, Dust, and Star-forming Properties (R1 - R39) 


ID 

Mh2 

^Paa 

SFR 

Tgas 

‘S '340 

-Mduat /M^ust 

Mism 


[xlO’^ Mq] 

[xlO^® erg s“l] 

[Mo yr-l] 

[Gyr] 

[mjy] 

[XIO'^ Mq] 

[xio’^ Mq] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) (8) 

(9) 


R1 

3.0 

± 

0.6 

< 24.6 


< 

0.15 


> 

0.20 


< 

0.51 


< 

17 


> 

180 


< 

3.6 

R2 

2.0 

± 

0.4 

46.0 

± 

8.2 

0.29 

± 

0.05 

0.07 

± 

0.02 


< 

0.51 


< 

17 


> 

119 


< 

3.6 

R3 

5.1 

± 

0.8 

11.0 

± 

8.2 

0.07 

± 

0.05 

0.75 

± 

0.57 


< 

0.51 


< 

17 


> 

306 


< 

3.6 

R4 

5.0 

± 

0.7 

9.3 

± 

8.2 

0.06 

± 

0.05 

0.88 

± 

0.78 


< 

0.51 


< 

17 


> 

299 


< 

3.6 

R5 

4.9 

± 

0.7 

41.8 

± 

8.2 

0.26 

± 

0.05 

0.19 

± 

0.05 


< 

0.51 


< 

17 


> 

293 


< 

3.6 

R6 

5.2 

± 

0.9 

60.7 

± 

8.2 

0.38 

± 

0.05 

0.14 

± 

0.03 


< 

0.51 


< 

17 


> 

310 


< 

3.6 

R7 

7.5 

± 

0.9 

172.1 

± 

8.2 

1.07 

± 

0.05 

0.07 

± 

0.01 

0.61 

± 

0.18 

20 

± 

6 

371 

± 

118 

5.2 

± 

1.4 

R8 

5.3 

± 

0.8 

69.8 

± 

8.2 

0.43 

± 

0.05 

0.12 

± 

0.02 


< 

0.51 


< 

17 


> 

314 


< 

3.6 

R9 

13.6 

± 

1.2 

125.1 

± 

8.2 

0.78 

± 

0.05 

0.18 

± 

0.02 


< 

0.51 


< 

17 


> 

811 


< 

3.6 

RIO 

18.0 

± 

1.3 

49.0 

± 

8.2 

0.30 

± 

0.05 

0.59 

± 

0.11 


< 

0.51 


< 

17 


> 1074 


< 

3.6 

Rll 

23.1 

± 

1.6 

106.9 

± 

8.2 

0.66 

± 

0.05 

0.35 

± 

0.04 


< 

0.51 


< 

17 


> 

1376 


< 

3.6 

R12 

6.5 

± 

0.9 

25.5 

± 

8.2 

0.16 

± 

0.05 

0.41 

± 

0.14 


< 

0.51 


< 

17 


> 

390 


< 

3.6 

R13 

7.4 

± 

1.0 

25.9 

± 

8.2 

0.16 

± 

0.05 

0.46 

± 

0.16 


< 

0.51 


< 

17 


> 

441 


< 

3.6 

R14 

6.7 

± 

0.8 

25.5 

± 

8.2 

0.16 

± 

0.05 

0.42 

± 

0.14 


< 

0.51 


< 

17 


> 

398 


< 

3.6 

R15 

7.1 

± 

0.9 

14.7 

± 

8.2 

0.09 

± 

0.05 

0.78 

± 

0.44 


< 

0.51 


< 

17 


> 

425 


< 

3.6 

R16 

7.9 

± 

0.9 

15.9 

± 

8.2 

0.10 

± 

0.05 

0.80 

± 

0.42 


< 

0.51 


< 

17 


> 

473 


< 

3.6 

R17 

9.1 

± 

1.0 

35.2 

± 

8.2 

0.22 

± 

0.05 

0.42 

± 

0.11 


< 

0.51 


< 

17 


> 

545 


< 

3.6 

R18t 

23.5 

± 

1.6 

277.4 

± 

8.2 

1.72 

± 

0.05 

0.14 

± 

0.01 

5.57 

± 

0.58 

183 

± 

19 

128 

± 

16 

47.1 

± 

4.9 

R19 

15.7 

± 

1.3 

150.8 

± 

8.2 

0.94 

± 

0.05 

0.17 

± 

0.02 

2.70 

± 

0.32 

89 

± 

11 

177 

± 

25 

22.9 

± 

2.6 

R20 

26.0 

± 

1.7 

78.1 

± 

8.2 

0.48 

± 

0.05 

0.54 

± 

0.07 

2.40 

± 

0.29 

79 

± 

10 

329 

± 

46 

20.3 

± 

2.4 

R2ltt 

46.9 

± 

2.7 

303.2 

± 

8.2 

1.88 

± 

0.05 

0.25 

± 

0.02 

8.55 

± 

0.87 

281 

± 

29 

167 

± 

19 

72.3 

± 

7.3 

R22 

40.7 

± 

2.4 

197.1 

± 

8.2 

1.22 

± 

0.05 

0.33 

± 

0.02 

4.08 

± 

0.44 

134 

± 

15 

303 

± 

37 

34.5 

± 

3.7 

R23 

47.8 

± 

2.7 

217.5 

± 

8.2 

1.35 

± 

0.05 

0.35 

± 

0.02 

3.60 

± 

0.40 

119 

± 

13 

403 

± 

50 

30.5 

± 

3.3 

R24 

22.9 

± 

1.7 

145.5 

± 

8.2 

0.90 

± 

0.05 

0.25 

± 

0.02 

1.28 

± 

0.21 

42 

± 

7 

543 

± 

99 

10.8 

± 

1.7 

R25 

35.8 

± 

2.2 

281.2 

± 

8.2 

1.74 

± 

0.05 

0.21 

± 

0.01 

3.10 

± 

0.35 

102 

± 

12 

351 

± 

46 

26.2 

± 

2.9 

R26 

29.0 

± 

1.9 

250.1 

± 

8.2 

1.55 

± 

0.05 

0.19 

± 

0.01 

1.66 

± 

0.24 

54 

± 

8 

533 

± 

84 

14.0 

± 

1.9 

R27 

17.7 

± 

1.4 

507.8 

± 

8.2 

3.15 

± 

0.05 

0.06 

± 

0.01 

1.42 

± 

0.22 

47 

± 

7 

379 

± 

66 

12.0 

± 

1.7 

R28 

8.2 

± 

1.0 

246.3 

± 

8.2 

1.53 

± 

0.05 

0.05 

± 

0.01 

0.73 

± 

0.18 

24 

± 

6 

344 

± 

97 

6.1 

± 

1.4 

R29 


< 

1.4 

147.8 

± 

8.2 

0.92 

± 

0.05 


< 

0.02 

0.76 

± 

0.19 

25 

± 

6 


< 

55 

6.4 

± 

1.4 

R30 

9.1 

± 

1.1 

40.2 

± 

8.2 

0.25 

± 

0.05 

0.36 

± 

0.09 

0.79 

± 

0.19 

26 

± 

6 

348 

± 

92 

6.7 

± 

1.4 

R31 

7.0 

± 

1.0 

19.0 

± 

8.2 

0.12 

± 

0.05 

0.59 

± 

0.27 


< 

0.51 


< 

17 


> 

418 


< 

3.6 

R32 

3.3 

± 

0.5 

60.3 

± 

8.2 

0.37 

± 

0.05 

0.09 

± 

0.02 


< 

0.51 


< 

17 


> 

198 


< 

3.6 

R33 

2.0 

± 

0.4 

55.0 

± 

8.2 

0.34 

± 

0.05 

0.06 

± 

0.02 


< 

0.51 


< 

17 


> 

122 


< 

3.6 

R34 

4.0 

± 

0.7 

10.9 

± 

8.2 

0.07 

± 

0.05 

0.59 

± 

0.45 


< 

0.51 


< 

17 


> 

238 


< 

3.6 

R35 

7.6 

± 

0.9 

< 24.6 


< 

0.15 


> 

0.51 


< 

0.51 


< 

17 


> 

453 


< 

3.6 

R36 

7.9 

± 

1.1 

< 24.6 


< 

0.15 


> 

0.52 


< 

0.51 


< 

17 


> 

469 


< 

3.6 

R37 

6.9 

± 

1.0 

< 24.6 


< 

0.15 


> 

0.46 


< 

0.51 


< 

17 


> 

410 


< 

3.6 

R38 

3.8 

± 

0.7 

15.2 ± 8.2 

0.10 

± 

0.05 

0.40 

± 

0.22 


< 

0.51 


< 

17 


> 

225 


< 

3.6 

R39ttt 

39.1 

± 

2.3 

273.6 ± 8.2 

1.70 

± 

0.05 

0.23 

± 

0.02 

3.50 

± 

0.58 

115 

± 

19 

339 

± 

60 

29.6 

± 

3.2 


Note. — Column 2: The molecular gas mass derived using the conversion factor aco = 0.8 (K km s“^ pc^)“^. Column 3: The 
Paa flux (Tateuchi et al. 2012). Column 4: The star formation rate derived using the conversion factor SFR/Lpaa = 6.2 X 10“^^ 
[erg s~^/(Mq yr~^)“^]. Column 5: The gas depletion time (= Sh 2 /^SPr)- Column 6: The 340 GHz continuum flux. Column 7: 
The dust mass derived using the equation (7). We adopt the Draine & Lee (1984) dust model for k; 34 o to derive the M^ust- Column 9: 
The ISM mass derived using the equation (8).; We consider the statistical error and the systematic error in this table. The systematic 
error of absolute flux calibration is estimated to be ~ 5 % in band 3 and ~ 10% in band 7. frepresents boxes contained the obscured 
AGN deflned by paper 1. tt^epresents boxes contained the nuclear starbursts defined by paper I. ftfrepresents boxes at the overlap 
region. 
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TABLE 11 

Gas, Dust, and Star-forming Properties (SO - S14) 


ID ^ dense ^Paa SFR, Tgas *^340 H^dust H^dense /-^dust 

[xlO® Mq] [xlO^S erg s-1] [Mq yr"!] [Myr] [mjy] [xlO^ M©] 

(1) (2) (3) (4) (5) (6) (7) (8) 


sot 

38.3 

± 

4.7 

322 

± 

12 

2.00 

± 

0.07 

19.2 

± 

2.5 

1.71 

± 

0.18 

56.1 

± 

6.0 

68 

± 

11 

SI 

3.2 

± 

1.8 

240 

± 

12 

1.49 

± 

0.07 

2.2 

± 

1.2 

0.25 

± 

0.07 

8.1 

± 

2.3 

40 

± 

25 

S2 


< 

4.3 


< 

36 


< 

0.21 




0.10 

± 

0.07 

3.1 

± 

2.2 

< 137 

S3tt 

13.9 

± 

2.6 

663 

± 

12 

4.11 

± 

0.07 

3.4 

± 

0.6 

1.28 

± 

0.14 

42.1 

± 

4.7 

33 

± 

7 

S4 

3.6 

± 

2.2 

446 

± 

12 

2.76 

± 

0.07 

1.3 

± 

0.8 

0.63 

± 

0.09 

20.6 

± 

3.0 

17 

± 

11 

S5tt 

23.2 

± 

3.6 

606 

± 

12 

3.75 

± 

0.07 

6.2 

± 

0.7 

2.56 

± 

0.26 

84.4 

± 

8.7 

27 

± 

5 

S6 

11.5 

± 

2.6 

268 

± 

12 

1.66 

± 

0.07 

6.9 

± 

1.6 

0.90 

± 

0.11 

29.7 

± 

3.7 

39 

± 

10 

S7 

5.4 

± 

2.3 

138 

± 

12 

0.85 

± 

0.07 

6.4 

± 

2.7 

0.40 

± 

0.08 

13.0 

± 

2.5 

42 

± 

19 

S8 

1.8 

± 

1.6 

158 

± 

12 

0.98 

± 

0.07 

1.8 

± 

1.6 

0.32 

± 

0.07 

10.5 

± 

2.4 

17 

± 

16 

S9 

2.0 

± 

1.8 

106 

± 

12 

0.66 

± 

0.07 

3.1 

± 

2.8 

0.28 

± 

0.07 

9.2 

± 

2.3 

22 

± 

21 

SIO 


< 

4.3 

203 

± 

12 

1.26 

± 

0.07 


< 

3.4 

0.34 

± 

0.07 

11.3 

± 

2.4 


< 

38 

Sllttt 


< 

4.3 

242 

± 

12 

1.50 

± 

0.07 


< 

2.9 

0.21 

± 

0.07 

7.0 

± 

2.3 


< 

61 

S12ttt 


< 

1.4 

323 

± 

12 

2.00 

± 

0.07 


< 

2.1 

0.46 

± 

0.08 

15.3 

± 

2.6 


< 

28 

S13 


< 

4.3 

291 

± 

12 

1.81 

± 

0.07 


< 

2.4 

0.46 

± 

0.08 

15.2 

± 

2.6 


< 

28 

S14 


< 

1.4 

219 

± 

12 

1.36 

± 

0.07 


< 

3.2 

0.27 

± 

0.07 

9.0 

± 

2.3 


< 

48 


Note. — Column 2: The molecular gas mass derived using the conversion factor ohcn — 10/0.63 (K km s ^ pc^) ^. 
Column 3: The Paa flux (Tateuchi et al. 2012). Column 4: The star formation rate derived using the conversion factor 
SFR/Lpaa — 6.2 X [erg s“^/(Mq yr“^)“^]. Column 5: The gas depletion time (— Sdense/^SFR)- Column 6 : 

The 340 GHz continuum flux. Column 7: The dust mass derived using the equation ( 8 ). We adopt the Draine Sz Lee 
(1984) dust model for K 340 to derive the Mdnst-'i We consider the statistical error and the systematic error in this table. 
The systematic error of absolute flux calibration is estimated to be ~ 5% in band 3 and ~ 10% in band 7. frepresents 
boxes contained the obscured AGN deflned by paper I. ffrepresents boxes contained the nuclear starbursts deflned by 
paper I. ftfrepresents boxes at the overlap region. 
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TABLE 12 

RADEX Results of Case 1 


ID 

(1) 

Tkin 
< 7.81 
[K] 

(2) 

min. x^ 

[K] 

(3) 

log 

X^ < 7.81 

(4) 

riH, 

min. X 

(5) 

loglV(H 2 ) 

X^ < 7.81 min. x^ 

[cm-2] [cm-2] 

(6) (7) 

min. x^ 

(8) 

R5 

> 5 

90 

< 3.7 

2.4 

20.4 - 21.2 

20.7 

0.0001508 

R7b 

5 - 120 

25 

> 2.8 

4.1 

> 20.8 

21.4 

0.0008348 

R8b 

5 - 30 

25 

< 2.4 

2.0 

20.9 - 21.8 

21.3 

0.04646 

R9b 

> 5 

50 

2.2 - 4.0 

3.0 

20.1 - 20.7 

20.4 

0.004822 

RlOa 

5 - 120 

65 

< 3.3 

2.1 

20.4 - 21.8 

20.9 

0.001529 

RlOb 

5 - 20 

10 

3.6 - 4.6 

4.4 

19.9 - 20.4 

20.1 

0.00948 

Rlla 

> 20 

75 

< 3.2 

2.5 

20.5 - 21.0 

20.8 

0.000315 

Rllb 


15 

2.6 - 4.6 

3.9 

< 20.4 

18.1 

0.001707 

R13a 

5 - 200 

85 

< 3.5 

2.2 

20.3 - 21.4 

20.7 

0.002488 

R14 

5 - 200 

15 

< 3.4 

3.2 

20.2 - 20.6 

20.4 

0,004464 

R15 

5 - 140 

35 

< 3.4 

2.6 

20.2 - 20.8 

20.5 

0.001343 

R16 

5 - 30 

15 

2.4 - 3.2 

2.8 

20.6 - 21.0 

20.8 

0.03139 

R17 

5 - 180 

135 

< 3.6 

2.0 

19.6 - 21.2 

20.2 

0.003854 

R18t 

> 35 

180 

2.4 - 3.3 

2.7 

20.6 - 21.0 

20.8 

0.000186 

R19a 

> 25 

85 

< 3.2 

2.5 

20.9 - 21.5 

21.2 

0.0001408 

R19b 

> 20 

100 

2.2 - 3.5 

2.8 

20.5 - 21.0 

20.8 

0.002962 

R20a 

> 30 

100 

2.7 - 3.6 

3.1 

20.5 - 21.1 

20.8 

0.001267 

R21att 

25 - 90 

50 

3.4 - 5.0 

3.7 

20.8 - 21.5 

21.1 

0.01197 

R22a 

> 25 

70 

2.6 - 3.6 

3.1 

20.5 - 20.9 

20.7 

0.01047 

R22b 

> 90 

265 

2.8 - 3.3 

3.0 

21.0 - 21.6 

21.3 

0.00147 

R23a 

> 40 

140 

2.7 - 3.5 

3.1 

20.6 - 21.1 

20.9 

0.001661 

R23b 

> 35 

140 

2.5 - 3.5 

2.9 

20.5 - 21.0 

20.8 

0.002198 

R24a 

> 5 

160 

< 4.9 

2.4 

20.2 - 21.3 

20.6 

0.000268 

R24b 

> 15 

275 

2.5 - 3.6 

2.8 

20.6 - 21.4 

21.0 

0.001229 

R24c 

> 5 

115 

< 3.9 

2.5 

20.0 - 20.8 

20.4 

0.0001756 

R25a 

5 - 40 

20 

3.5 - 5.0 

3.9 

20.5 - 20.8 

20.7 

0.002611 

R25b 

> 5 

285 

< 3.9 

2.2 

20.2 - 21.0 

20.4 

0.0001641 

R26a 

> 40 

265 

2.6 - 3.4 

2.8 

20.5 - 21.1 

20.8 

0.0002885 

R26b 

10 - 100 

50 

< 3.3 

2.2 

20.6 - 21.4 

20.9 

0.00003488 

R27a 

> 5 

255 

> 2.9 

3.0 

< 21.6 

20.6 

0.0004384 

R27b 

> 50 

160 

< 2.6 

2.0 

21.0 - 21.7 

21.3 

0.001546 

R28a 

> 5 

110 

< 3.9 

2.5 

20.6 21.6 

21.0 

0.00007702 

R31a 

> 20 

115 

< 3.0 

2.0 

20.8 - 22.0 

21.3 

0.001847 

R34 

> 5 

130 

< 3.4 

2.4 

20.5 - 21.6 

20.9 

0.0006237 

R36a 

10 - 40 

30 

< 2.5 

2.0 

> 21.1 

21.6 

0.00002201 

R37a 

5 - 30 

25 

< 2.7 

2.0 

20.7 - 21.8 

21.1 

0.03561 

R38 

> 5 

35 

< 3.7 

2.5 

> 20.8 

21.2 

0.0005419 

R39attt 

> 50 

95 

2.7 - 3.4 

3.1 

20.9 - 21.3 

21.1 

0.006462 

R39bttt 

> 5 

175 

2.3 - 4.1 

2.6 

20.0 - 20.9 

20.5 

0.005283 


Note. — Column 2-7: RADEX parameters noted above. Column 8: Value of associated with the 
fit. Note that for each position, (3), (5), and (7) are the best-fit parameters and (2), (4), and (6) are 
estimated within a confidence of 95 %.; We consider the statistical error and the systematic error in this 
table. The systematic error of absolute flux calibration is estimated to be ~ 5% in band 3 and ~ 10% in 
band 7. frepresents boxes contained the obscured AGN defined by paper I. ffrepresents boxes contained 
the nuclear starbursts defined by paper I. ftfrepresents boxes at the overlap region. 


TABLE 13 

RADEX Results of Case 2 


ID 

(1) 

logAr(H2) 

[cm-2] 

(2) 

Tkin 

x2 < 7.81 
[K] 

(3) 

min. x^ 

[K] 

(4) 

log 

X^ < 7.81 
[cm-3] 
(5) 

riHa 

min. x^ 

[cm-3] 

(6) 

[HCN]/[HCO+] 

X^ < 0.35 min. x^ 

[cm-2] 

(7) (8) 

min. x^ 

(9) 

EOt 

21.2 

> 100 

270 

5.0 - 5.4 

5.3 

> 5 

8 

0.00563 

Eltt 

21.6 

40 - 100 

70 

5.6 - 5.9 

5.8 

< 4 

2 

0.01073 

E2ttt 

21.5 

5 - 90 

40 

5.0 - 5.6 

5.1 

1 - 9 

2 

0.0003303 


Note. — Column 2: Adopted A(H 2 ) which derived from CO (1-0) data and Xco- Column 3-8: RADEX 
parameters noted above. Column 9: Value of associated with the fit. Note that for each position, (4), (6), 
and (8) are the best-fit parameters and (3), (5), and (7) are estimated within the confidence level of 95 % for Tkin 
and 71^2 y ci-nd 95 % (2 a) for [HCN]/[HCO'*’].; We consider the statistical error and the systematic error in this 
table. The systematic error of absolute flux calibration is estimated to be ~ 5% in band 3 and ~ 10% in band 
7. frepresents boxes contained the obscured AGN defined by paper I. ffrepresents boxes contained the nuclear 
starbursts defined by paper I. fffrepresents boxes at the overlap region. 
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TABLE 14 

Gas Properties under the LTE Assumption 


ID 

(1) 

'S’A (— Tb_cO(l-0)/Tkin) 

(2) 

Adopted $A 

(3) 

Tex 

[K] 

(4) 

■^00(1-0) 

(5) 

■^1300(1-0) 

(6) 

R5 

< 0.24 

0.10 

15.3 ± 0.9 

3.59 

0.14 

R7b 

< 0.09 

0.09 

8.3 ± 0.8 

3.97 

0.02 

R8b 

0.03 - 0.19 

0.10 

13.0 ± 0.9 

25.28 

1.57 

R9b 

< 0.52 

0.10 

29.5 ± 1.5 

1.27 

0.02 

RlOa 

0.01 - 0.16 

0.10 

11.4 ± 0.8 

6.72 

0.35 

RlOb 

0.14 - 0.56 

0.14 

23.5 ± 1.1 

1.24 

0.03 

Rlla 

< 0.11 

0.10 

24.6 ± 1.3 

3.88 

0.19 

Rllb 


0.10 

36.9 ± 1.8 

0.01 

0.01 

R13a 

0.01 - 0.22 

0.10 

14.6 ± 0.9 

6.87 

0.37 

R14 

0.01 - 0.32 

0.10 

19.6 ± 1.1 

5.89 

0.35 

R15 

0.01 - 0.35 

0.10 

21.2 ± 1.1 

6.44 

0.27 

R16 

0.05 - 0.28 

0.10 

17.2 ± 1.0 

9.00 

0.31 

R17 

0.01 - 0.33 

0.10 

19.9 ± 1.1 

2.59 

0.06 

R18t 

< 0.10 

0.10 

38.7 ± 1.9 

1.34 

0.03 

R19a 

< 0.04 

0.04 

28.7 ± 2.2 

6.99 

0.76 

R19b 

< 0.10 

0.10 

23.4 ± 1.2 

1.68 

0.05 

R20a 

< 0.11 

0.11 

36.8 ± 1.8 

0.57 

0.01 

R21att 

0.06 - 0.22 

0.10 

58.8 ± 2.9 

2.75 

0.03 

R22a 

< 0.14 

0.10 

39.5 ± 1.9 

0.90 

0.01 

R22b 

< 0.04 

0.04 

99.3 ± 5.1 

0.42 

0.01 

R23a 

< 0.10 

0.10 

45.2 ± 2.2 

0.61 

0.01 

R23b 

< 0.10 

0.10 

37.4 ± 1.8 

0.87 

0.01 

R24a 

< 0.25 

0.10 

15.7 ± 0.9 

4.00 

0.21 

R24b 

< 0.11 

0.10 

19.5 ± 1.1 

1.02 

0.01 

R24c 

< 0.36 

0.10 

21.3 ± 1.1 

1.83 

0.05 

R25a 

0.12 - 0.98 

0.12 

44.3 ± 2.1 

2.85 

0.03 

R25b 

< 0.35 

0.10 

21.2 ± 1.1 

2.16 

0.07 

R26a 

< 0.11 

0.10 

48.4 ± 2.4 

1.56 

0.02 

R26b 

0.01 - 0.11 

0.10 

14.1 ± 0.9 

6.33 

0.30 

R27a 

< 0.39 

0.10 

23.1 ± 1.2 

0.09 

0.01 

R27b 

< 0.03 

0.03 

54.5 ± 3.5 

9.17 

0.92 

R28a 

< 0.15 

0.10 

11.0 ± 0.8 

4.21 

0.31 

R31a 

< 0.04 

0.04 

21.5 ± 2.0 

4.61 

0.09 

R34 

< 0.17 

0.10 

11.7 ± 0.8 

4.23 

0.29 

R36a 

0.02 - 0.07 

0.07 

12.7 ± 1.1 

20.17 

1.53 

R37a 

0.03 - 0.16 

0.10 

11.3 ± 0.8 

19.92 

1.04 

R38 

< 0.08 

0.08 

8.0 ± 0.9 

10.07 

0.67 

R39attt 

< 0.10 

0.10 

52.6 ± 2.6 

1.42 

0.01 

R39bttt 

< 0.39 

0.10 

23.1 ± 1.2 

1.29 

0.02 


Note. — Column 2: The beam filling factor estimated from RADEX modeling and the brightness 
temperature. Column 3: The adopted beam filling factor. Column 4: The excitation temperature of 
a spectral line calculated from the equation (3). Column 5 and 6: the optical depth of ^^CO (1—0) and 
^^CO (1—0) emission from RADEX modeling, respectively. Column 7: The total column density of 
^^CO.; We consider the statistical error and the systematic error in this table. The systematic error 
of absolute flux calibration is estimated to be ~ 5% in band 3 and ~ 10% in band 7. frepresents boxes 
contained the obscured AGN defined by paper I. ffrepresents boxes contained the nuclear starbursts 
defined by paper 1. ftfrepresents boxes at the overlap region. 
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TABLE 15 

Properties of CS, CH 3 OH, and CN 


Molecule 

( 1 ) 

ID 

( 2 ) 



Tb 

[K] 

(3) 

At; 
[km s“^] 

(4) 



Nx 

[cm 

(5) 

{x]/m 

( 6 ) 

CS 

RlSt 


< 

0.04 


< 

4.0 

X 

lO^i 

< 6.3 X 10-1° 


R21att 


< 

0.04 


< 

5.6 

X 

loll 

< 4.5 X 10-1° 


R39attt 

0.13 

± 

0.04 

44 ± 11 


1.6 

X 

1012 

1.3 X 10-° 

CH 30 H 

RlSt 


< 

0.04 


< 

1.9 

X 

1012 

< 2.9 X 10-° 


R21att 


< 

0.04 


< 

2.6 

X 

1012 

< 2.1 X 10-° 


R39attt 

0.17 

± 

0.04 

43 ± 5 


1.0 

X 

1013 

8.1 X 10-° 

CN 

RlSt 

0.23 

± 

0.04 

63 ± 5 


2.4 

X 

1012 

3.8 X 10-° 


R21att 

0.15 

± 

0.04 

66 ± 3 


2.2 

X 

1012 

1.8 X 10-° 


R39attt 

0.14 

± 

0.04 

49 ± 9 


1.4 

X 

1012 

1.1 X 10-° 


Note. — Column 1: The molecular line considered. Column 3: Peak brightness 
temperature of considered line. Column 4: Line width of considered line. Column 5: 
the derived total box-averaged column density of this species. Column 6: Fractional 
abundance relative to H 2 of this species.; We consider the statistical error and the 
systematic error in this table. The systematic error of absolute flux calibration is 
estimated to be ~ 5% in band 3 and ~ 10% in band 7. frepresents boxes contained the 
obscured AGN defined by paper I. ft^’^P^^sents boxes contained the nuclear starbursts 
defined by paper 1. ftfrepresents boxes at the overlap region. 
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APPENDIX 

A. IMAGES AND SPECTRA OF VV114 
A.l. The Paa and Kg band image with miniTAO/ANIR observation 



Fig. 14.— The Pao and Kg band images of VV 114 with miniTAO/ANIR observation (Tateuchi et al. 2012). 
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A.2. Channel maps of each line emission 



Fig. 15.— The velocity channel maps of the (1—0) line emission of VV 114. Two crosses in each channel show the positions of the 

nuclei defined by the peak positions of the Ks-band observation (Tateuchi et al. 2012). The velocity width of each channel is 10 km s“^. 
The beam size is plotted in the bottom-left corner of each channel. The contours represent flux intensity levels: -4.6, 4.6, 9.2, 18.4, 36.8, 
73.6, 110.4 and 147.2 mjy beam”!. 
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Fig. 15.— Continuied. 
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Fig. 16.— The same as Figure 14 but for the CO (3-2) line emission of VV 114. The approximate field of view of ALMA 7-point mosaic 
at this frequency is indicated by the large red circle. The contours represent fiux intensity levels: -12.6, 12.6, 25.2, 50.4, 100.8, 201.6, and 
403.2 mjy beam~^. 
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Fig. 16.— Continuied. 
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Fig. 17.— The same as Figure 14 but for the HCN (4-3) line emission of VV 114. The velocity width of each channel is 30 km s ^. The 
cross in each channel shows the position of the eastern nucleus defined by the peak position of the Ks-band observation (Tateuchi et al. 
2012). The contours represent flux intensity levels: -2.4, 2.4, 4.8, 9.6, and 19.2 mjy bearn”^. 
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Fig. 18.— The same as Figure 14 but for the HCO"*' (4-3) line emission of VV 114. The velocity width of each channel is 30 km s“^. 
The cross in each channel shows the position of the eastern nucleus defined by the peak position of the Ks-band observation (Tateuchi et 
al. 2012). The contours represent flux intensity levels: -2.4, 2.4, 4.8, 9.6, and 19.2 mjy bearn”^. 
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A.3. Box-summed spectra of each line emission 
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Fig. 19.— 2'.'0 box-summed spectra of (1-0) (red line), (1-0) X 10 (green broken line), and (3-2) (blue broken line) 

at the each box, labeled R1 - R39 of Fig 7. The spectra are taken from the ALMA data cubes after correcting the cubes for the primary 
beam attenuation and convolving them to 2'.'0 X I'.'S resolution (P.A. = 83 deg). 
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Fig. 19.— continued. 



Fig. 19.— continued. 
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Fig. 20.— 2^0 box-summed spectra of (3-2) (red line), CN (l 3 / 2 “ 0 i/ 2 ) 10 (green broken line), CS (2-1) x 10 (blue broken line), 

and CH 3 OH (2fc-lfc) x 10 (light blue broken line) at R18, R21, and R39, of Fig 7. The spectra are taken from the ALMA data cubes after 
correcting the cubes for the primary beam attenuation and convolving them to 270 x l75 resolution (P.A. = 83 deg). 


Fig. 21.— 
Fig 7. 



C/.'5 box-summed spectra of HCN (4-3) (red line) and HCO"*” (4-3) (green broken line) at the each box, labeled SO - S14, of 





















































































































































